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ABSTRACT 



We have developed a self-co n sistent - field technique similar to the one described 
bv lHachisu. Eriguchi. fc Nomotd (jl986bl ) that can be used to construct detailed force- 
balanced models of synchronously rotating, double white dwarf (DWD) binaries that 
have a wide range of total masses, mass ratios, and separations. In addition to providing 
a computational tool that can be used to provide quiet initial starts for dynamical studies 
of the onset of mass transfer in DWD systems, we show that this SCF technique can be 
used to construct model sequences that mimic the last portion of the detached inspiral 
phase of DWD binary evolutions, and semi-detached model sequences that mimic a 
phase of conservative mass transfer. 



Subject headings: binaries: close (AM CVn), hydrodynamics, methods: numerical, 
white dwarfs, supernovae 



1. Introduction 



As IWebbinkI (j 19841 ) and llben fc Tutukovl (|l984 119861 ) have pointed out, double white dwarf 
(DWD) binaries are expected to be the end product of the thermonuclear evolution of a sizeable 
fraction of all binary systems. The subset of DWD binaries that are born with orbital periods 
-forb ^ 5 hours are of particular astrophysical interest because they will be driven in to contact 
withi n a Hubble time via the loss of angular momentum through gravitational radiation (jPacvzhski 
19671 ). Even this short-period component of the DWD popula tion is expected to be quite lar ge in our 
Galaxy. Sufficiently large that, according to lHils et al.l (| 19901 ) and lCornish &: LarsonI (|2003l ). DWDs 
are likely to be a dominant source of background noise for t he proposed .Laser Interferometer Space 
Antenna [LISA) (jFaller &: Bendeiill984l : lEvans et al.lll987l : lBenderlll998l ) in its lower gravitational- 
wave frequency band, /gw = 2/Porb ^ 3 x 10~^ Hz. Ironically, it is difficult to detect DWDs via 
traditional observational techniques because the intrinsic photon luminosity of white dwarfs is very 
low. 

Broadly speaking, our Galaxy's DWD binary population should be dominated by systems that 
are in two distinctly different evolutionary phases: the inspiral phase alluded to above, during which 
both stars are detached from their respective Roche lobes; and a semi-detached, stable mass transfer 
phase, during which the less massive star fills its Roche lobe and is slowly trans ferring mass to its 
more massive companion. To date ~ 100 detached systems have been identified (iNapiwotzki et al 
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200ll . |2004| ). Orbital pe riods and component masses have been determined for approximately one 
quarter of this samp le (|Marshl [2OO0I : iMaxted et aPboOol : iNapiwotzki et alibooa : Ixarl et alJbood : 
Nelemans et al.ll2005l ) . AM CVn is the prototype of semi-detached DWD systems that are undergo- 
ing a phase of stable mass transfer (IWarneiill995l). In the immediate solar neighborhood ^ 20 suc h 



2005 



Ramsay et ahlboO?! ). 



syste r ns are known (IWoudt k. Warneij l2003l : iNelemana l2005l : lAnderson et al.l l2005l : iRoelofs et al 



Theoretical arguments suggest that the ultimate fate of a DWD binary will depend on the 



system's total mass, Mtpt, a r id mass rat io go at the onset of mass transfer (jMarsh et al.l l2004l : 



Gokhale. Peng &: Frank! 120071 : iFrankI 120081 ) . For example, a DWD will likely only be able to enter 
an extended phase of stable mass transfer as characterized by AM CVn systems if is less than 
some critical value ^stable — where (?stabie is almost certainly < 2/3 and may be closer to 1/4. If 
^0 > (Zstabie! the mass transfer rate is expected to diverge in a finite time, ultimately implying tidal 
disruption of the donor and/or a catastrophic rn erger of the two components. Consiste nt with 
the theoretical ideas presented by IWebbinkI (119841 ) and llben. Tutukov &: YungelsonI (jl996l ). recent 
observations strongly suggest that the end product of so me DWD mergers a re R Coronae Borealis 
(RCB) stars and hydrogen-deficient carbon (HdC) stars (jClavton et al.ll2007l ). Also, DWD systems 
with Mtot in e xcess the Chandrasekhar mass have long been considered like l y progenitors of T ype 



la super novae (jWebbinklll984l : llben. Tutukov fc Yungelson 



1996 



Livio 



2000 



YoonetalJ 12003). 



Over the past couple of decades, various groups have employed smoothed particle hydrody- 

entl 



namics (SPH) techniques to illustrate the dynamical be havior of DWD systems that viq ^ 
merge after encountering an unstable mass-transfer event ( Benz et al. 199ol: Rasio &: Shapirol 



Yoon et al 



Segretain. Chabrier. &: Mochkovitch 1997 : Fryer et al. 19991 : Guerrero. Garcia-Berro. fc IsernI 



entiy 
1993 : 



200 



1; 



20071 ). Typically, initial states for these simulations have been constructed in such a 



way that the rnerger process is completed in < 5 orbits after first contact. iD'Souza et al.l (|2006l ) 
and iMotl et al.l (j2007l ) recently employed a grid-based finite- volume (FV) computational fluid dy- 
namic technique to also study the onset and nonlinear development of mass transfer in strongly 
interacting binary systems. Their primary objective was to ascertain the value of ^stable in sys- 
tems, such as DWDs, that undergo a phase of direct impact accretion following the onset of mass 
transfer. They were able to follow the evolution of a couple of different systems through > 30 
orbits. Instead of merging, these systems appeared to be entering a long-term phase of stable mass 
transfer. Presumably these are the types of binary configurations that serve as the progenitors of 
AM CVn systems. 



Motl et al.l (|2007l ) found that the outcome of their simulations — for example, whether a system 
survives the onset of mass transfer or merges — can be sensitive to initial conditions. In particular 
if, rather than making only marginal contact with its Roche lobe, the donor star is in relatively 
deep contact at the onset of a simulation, the mass-transfer rate will initially be artificially high 
and it may be difficult for the system to avoid merger. As a result, a DWD system that should be 
categorized by numerical simulation as an AM CVn progenitor may be incorrectly categorized as 
a likely progenitor of an RCB star or a Type la supernova. It is therefore important to start such 
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simulations from initial states that are very quiet — that is, from initial configurations that are 
in detailed force balance throughout — and to perform each simulation with a grid (or particle) 
resolution that is sufficient to resolve margina l, or at least ve r y sha llow, contact between the donor 
and it s Roche lobe. In two very recent reports. iFrver &: Diehll (|2008l ) and lDan. Rosswog. &: Briiggen 
(|2008l ) have confirmed this finding. Both groups have shown that SPH techniques also can be used 
to follow stable mass-transfer events through > 30 orbits if sufficiently quiet initial states are used 
and if the simulations are carried out using a sufficiently large number of SPH particles. 



Quiet i nitial states were co nstructed for the hydrodynamic simulations reported bv lD'Souza et al 



(j2006l ) and iMotl et al.l (120071) us ing a self-consistent-field (SCF) technique very similar to the 
one described bv iHachisul (|l986bl ). However, these initial models were constructed using a poly- 
tropic equation of state. He r e we u se a modified version of the SCF technique developed by 
Hachisu. Eriguchi. &: Nomotd (1986bl) for the zero temperature white dwarf (ZTWD) equation of 



state derived by IChandrasekhaii (119351 ). We show how this SCF technique can be used to construct 
synchronously rotating, unequal-mass DWD binaries in which the less-massive (donor) star is in 
marginal contact with its Roche lobe and thereby provide excellent initial models for dynamical 
studies of interacting DWD systems that have a realistic equation of state. 

To demonstrate the broad utility of this SCF technique, we construct a sequence of detached, 
synchronously rotating binaries of varying separation, but fixed mass ratio and constant total mass, 
to mimic the last portion of the inspiral phase of evolution of DWD binary systems. This enables us 
to determine, for example, the degree to which the functional dependance of the orbital frequency 
on orbital separation 0(a) deviates from a pure Keplerian behavior. This type of sequence can be 
used to accurately identify the critical separation Ocrit and corresponding orbital frequency at which 
the less massive component of a DWD binary first makes contact with its Roche lobe. Information 
of this nature will be helpful in inte rpreting future observa t ions t hat are expected to be made with 
LISA, as described for example bv iKopparapu fc Tohlind (|2007l ). Finally, we show how this SCF 
technique can be used to construct a sequence of semi-detached binaries of fixed total mass, but 
varying separation and decreasing mass ratio, to mimic the initial stage of evolution of a DWD 
binary that has entered an AM CVn (i.e., stable mass-transfer) phase of evolution. 



2. SCF Formulation 



The SCF technique was first introduced to the astrophysics community by lOstriker &: Mark 



(| 19641 ) to create models of rapidly rotating, single stars with a polytropic equation of state. Hachisu 
developed a variation of the technique, improving convergence rates and extending its capabilities 
to include the use of a ZTWD equation of state. With his improved technique, Hachisu was 
able to construc t two-dimensional (2D) configurations of differentially rotating, single white dwarfs 
([Hachisul Il986al ) and three-dimensional (3D) config urations of unifo r mly rotating mu l tiple white 
dwarf systems in which the stars have equal mass (jHachisul Il986bl ). iNew fc Tohlind ()1997l ) em- 
ployed Hachisu's 3D technique to construct inspiral sequences of equal-mass DWD binaries, includ- 
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ing over-contact models having separations even smaller than Ocrit- iHachisul (|l986bl ) also applied 
his technique to the construction of unequal-mass binary systems using a polytropic equation of 
state and after additional algorithmic innovations were introduced, Hachisu, Eriguchi, and Nomoto 
(1986a, b) constructed a small sample of unequal-mass DWD binaries and heavy-disk white dwarf 
systems to examine the likely outcome of DWD mergers. In what follows we show how Hachisu's 
SCF technique for constructing unequal-mass DWD binaries can be further improved and used to 
construct inspiral binary sequences. 



2.1. Equation of State 



In the zer o-temperature white dwarf (ZTWD) equation of state (jChandrasekhar 



1935 



1967 



Hachisulll986al ) the electron degeneracy pressure P varies with the mass density p according to the 
relation, 



A 



where the dimensionless parameter, 



;(2x2 - 3)(x2 + 1)^/2 + 3sinh-i 



(1) 



1/3 



and the constants A and B are (see Appendix |A] and Table [AD for elaboration), 



4 5 

A = = 6.00228 X 10^^ dynes cm'^ , 



B 



3/l3 



9.81011 X 10^ g cm"^ . 



(2) 

(3) 
(4) 



According to IChandrasekhaij (|l967l ) (see again our Appendix [X]) , a natural length scale associated 
with models of ZTWDs is, 



/ 2A \ ^^"^ u 

pj^ = f—j ^ = 7.71395 X 10^ cm = O.Ollli?0 , 



(5) 



and the associated limiting white dwarf mass is 

3/2 / .. \ 2 



47r(2.01824) (j^^ (^^) = 1-14205 x lO^'' g = 5.742M0 



(6) 



Throughout this work, we will assume that the average ratio of nucleons to electrons throughout 
each white dwarf is = 2. Hence, B = 1.96202 x 10^ g cm" -3, = 5.55 X IQ-^Rq, and M^h = 
1.435Mo. 
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In terms of the enthalpy of the gas0 



H 



dP 
P 



the ZTWD equation of state shown in Eq. ([T|) can also be written in the form, 



H 



B 



+ 1 



1/2 



Inverting this gives the dependence of p on namely, 

„3 



P_ 
B 



- 1 



V 8^ 



(7) 



(8) 



(9) 



As a foundation for both constructing and understanding the structures of the synchronously 
rotating and tidally distorted stars in ZTWD binary systems, we have regenerated Chandraskehar 's 
spherical white dwarf sequence using a variation of the SCF technique outlined bv lHachisul (|l986al ). 
As is discussed in ^3.1\ Table [2] details key properties of the ZTWD structures that lie along this 
spherical model sequence. The white dwarf mass-radius relationship that is derived from models 
along this sequence is illustrat ed by the diamo nds in Figure [H For comparison, results from the 
published spherical sequence of iHachisul (jl986al ) are represented in this figure by asterisks and the 
solid cu rve shows the appr oximate, analytic mass-radius relationship, Eq. (jA14p . derived for ZTWD 
stars by lNauenberg] (119721 ). (As explained in Appendix El it is more appropriate for us to compare 
our results to this "Nauenberg" mass-radius relation than to the more widely used "Eggleton" 
mass-radius relation, shown in Eq. IA16[ ) 



2.2. Binary System Geometry and Governing Equations 

Our objective is to determine the 3D structure of a pair of ZTWD stars that are in a tight, 
circular orbit under the condition that both stars are synchronously rotating with the binary 
orbital frequency Q. We begin by specifying the masses Mi and M2 of the primary and secondary 
stars, respectively, such that M2 < Mi. Alternatively, we can specify the total system mass 
Mtot = Ml + M2 and the system mass ratio q = M2/M1 < 1, in which case, 

Ml = (^)m.„ 
M2 = (:r^lMtot 



Figure [2] shows a slice through the equatorial plane of such a system under the assumption that 
both stars are spherically symmetric. For our final equilibrium models in which the effects of 



^As defined here, H is actually enthalpy per unit mass. 
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tidal and rotational distortions are taken into account in a fully self-consistent fashion, this figure 
provides only a schematic illustration of the binary system's equatorial-plane structure. However, 
it provides an accurate depiction of the equatorial-plane structure of the initial stellar models that 
are fed into our iterative SCF scheme (see ^2.3p . 



In Figure [21 the more massive, primary star is shown on the left and the less massive, secondary 
star is on the right; the centers of the stars are located a distance w\ and respectively, from the 
center of mass of the binary system; and the binary separation a = vo\ + vo^- Because we are using 
a ZTWD equation of state, the central density /oJnax ^'^d radius Ri of each star (i = 1,2) cannot 
be specified independently of each star's chosen mass. As an initial guess for our SCF technique, 



the values of p 



max' rmax 



, and i?2 are drawn from Tabled that is, they are given by values that 



correspond to spherical ZTWDs having masses Mi and M2. 

For various values of the three principal system parameters Mi, M2 and a, our specific aim 
is to determine in a self-consistent fashion on a cylindrical coordinate mesh {w,9,Z), the values 
and spatial distribution of the scalar fields p^{zu, 6, Z), W[w, 6, Z), and P^{zu, 6, Z), for both stars 
(i = 1,2) in synchronously rotating, ZTWD binaries. Following Hachisu ( 1986a, b) and Hachisu, 
Eriguchi, and Nomoto ( 1986a, b), in order to construct these desired binary configurations we need 
to solve the following five equations simultaneously: 



p {w, 



47rG^p^(ti7, 



H\w,e,z) + ^{w,e,z) 



B[x\w,e,z)f = B 



-\ 2 



1 



3/2 



(10) 

(11) 
(12) 



where $(ro, 6, Z) is the Newtonian gravitational potential of the combined stellar system, and C* ^ 
and C*^^ are constants that specify the conditions of the Bernoulli flow inside each star. 



2.3. Solution Strategy 



Our solution strategy follows closely the methods described in lHachisul (jl986bl ) and lHachisu. Eriguchi. Nome 



(jl986bl ) so only the differences between our methods will be given in detail here. The variables with 
a carat (') abov e them are the dimensionless version of the variables as defined in Eqs.(22)-(27) of 
Hachisul (|l986bl ^. 



To begin an SCF iteration two spherical stars are initially placed on the computational grid 
in such a way that the center of mass of the system falls at the origin of the coordinate system and 
the outer edge of the secondary star (point Oa in Figure [2]) is at tz7 = 1. The centers of the stars 
are therefore located, respectively, at 

CTi _ q 



(13) 
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W2 = 



TUn 



1 



(14) 



where the dimensionless ratio i = R2/a is known once M2 and a have been chosen. These two 
expressions make sense because Wa = {vJ2 + o,tj = \w2 + {j^2 + and, for a point-mass binary 

whose center of mass is at the origin of the grid, vox = qw2- With defined everywhere on the 
grid, $(zz7,0,z) is calculated via Eq. ([TO]) . In this work the boundar y values for ^ ar e calcu lated 
using the compact cylindrical Green's function expansion described in lCohl &: Tohlind (|l999l ). and 
the values of the potential throughout the interior volume of the computation al grid are calcula ted 
using the Krylov subspace methods provided by the PETSc software library (jBalay et al.ll2004l ). 



During each iteration, t he interior structure of the secondary s tar is calculated using the same 
strategy as outlined in §2b of lHachisu. Eriguchi. &: Nomotd (jl986bl ). To update the structure of the 
primary star, h owever, the location of H'^^}^ is used in stead of specifying the inner edge of the star, 
as was done in iHachisu. Eriguchi. &: Nomotd (|l986bl ) . This choice eliminates the extra iterative 
steps that were needed when the third boundary condition was specified at the inner edge of the 
primary star. Using Eq. ([8]), we determine the value of the normalized enthalpy at the center of 
the primary star from the values of x^^^-^ selected for both stars and the value of H^^^ just derived 
for the secondary star. Specifically, we set 



HI 



1/2 



(15) 



Then we determine the value of the Bernoulli consta nt inside the pri mary star by examining the 
distribution of the variable F as defined in Eq.(48) of iHachisul (|l986al ). that is. 



4 + -(fw"^ 
2 



(16) 



In the vicinity of the original center of the primary star, that is, in the vicinity of point Oi as 
illustrated in Figure [2l the function F should exhibit a local maximum. We associate the location 
of this local maximum with the updated position of point Oi and we set -F^^ax equal to the value 
of the function at this local maximum. We therefore deduce from Eq. (jl6p that. 



HI 



F' 



(17) 



With this constant in hand, the normalized enthalpy throughout the primary star can be determined 
via the expression. 



1=1 



(18) 



and we obtain an updated "guess" for the normalized density distribution inside the primary star 
via the expression. 



sj=l 



1 



(xj^2 )3 1 \jji=l 

V max/ y. x-'-'ma,^ 



fri=l \ 2 >^ 3/2 

[«=ai)' + 1] - 1 



(19) 
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2.4. Global Properties and Convergence 



Our iterative scheme is judged to be operating well if various calculated model parameters 
— such as the dimensionless stellar masses M* and Bernoulli constants C* — converge toward 
well-defined values. We also have found it useful to track the convergence of various global energy 
parameters. Specifically, at the end of each iteration cycle we calculate the dimensionless rotational 
kinetic energy gravitat ional potential energy total internal energy U (see, for example, 
Eq. (75') in Chapter XI of I Chandr asekharl (jl967l )). and globally averaged pressure IT of the model, 
defined as follows: 



K 



W = 



U = 



n 



-<PpdV , 



H 



PdV . 



8A\ . 



dV, 



(20) 
(21) 
(22) 
(23) 



where dV = wdwdOdz is the dimensionless differential volume element on our cylindrical grid. 
Then the system's dimensionless total energy is given by the sum. 



Ei 



tot 



tot 



i=2 ^2^5 
max/ 



K + W + U. 



(24) 



and, if the model has converged to a proper equilibrium state, according to the virial theorem we 
should expect. 



2K + W + 3n 



0, 



(25) 



In general, at each iteration step the condition of virial equilibrium, Eq. (j25p . will not be satisfied, 
but if our iteration scheme is well behaved, convergence toward the virial condition should be 
achieved. With this in mind, we have found that the virial error. 



VE = 



2K + W + 3n 



W 



(26) 



provides a meaningful measure of the quality of each model. 



We declare that satisfactory convergence to a given model has been achieved when the abso- 
lute value of the fractional change between iterations has dropped below a specified convergence 
criterion, 5 ~ 10~^, for all of the following quantities: (7*, M*, Cl, K, W, 11, and the physical 
value of -cja. In addition, the converged model is judged to be a good equilibrium state if the virial 
error, VE, is sufficiently small. Table [J illustrates how we were able to achieve a lower virial error 
and, hence, a more accurate representation of an equilibrium configuration, by improving the grid 
resolution and/or by specifying a tighter convergence criterion. Specifically, the table shows that as 
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we were constructing binary model B3 (see discussion associated with Table [6l below) we were able 
to push the VE down from a value ~ 5 x 10~^ to a value ~ 5 x 10~^ by increasing the grid resolution 
from (64,128,33) to (128,256,65) zones in {wfi,z) and by pushing 5 from lO'^ to 3.5 x IQ-^. 

After the SCF code has converged to the desired equilibrium model, the various dimensionless 
variables are conver ted back t o prop er physical units following, for example, the scalings presented 
in Eqs.(22)-(27) of iHachisul (|l986al ). We note in particular that the value of the scale length 
vjt, = Wa is obtained by evalu ating a dimensio nless version of Eq. ([8|) for the secondary star in 
combination with Eq.(27) from IHachisul (|l986al ). which gives. 



8A/B 



,1=2 
max 



1/2 



-1/2 



(X 



i=2 \2 



+ 1 



1/4 



(27) 



In addition to the physical variables already identified, for each converged model we have found it 
useful to evaluate the system's total angular momentum. 



tot 



(28) 



as well as the spin angular momentum of each component star, Jg*pj„, and each star's Roche- 
lobe filling factor, f^^. As with the determination of quantities such as M* and R^, these latter 
two quantities are obtained by performing volume integrals over appropriate sub-domains of the 
computational grid, determined as follows. Let the origin of a Cartesian grid coincide with the 
center of mass of the binary system and align the x-axis of that grid with the line that connects 
the centers of the two stars as illustrated in Figure [2j Between points Oi and O2 along this axis, 
the effective potential. 



$cfr(x) = $(x) 



(29) 



will exhibit a maximum at position x^i associated with the inner "LI" Lagrange point. We define 
sub-domain P^^^ as the volume of the grid for which x = wcosO > xli and p > 0, that is, the 
region occupied by the secondary star; and we define sub-domain PJ,^^ as the volume of the grid 
for which x < x^i and p > 0, that is, the region occupied by the primary star. Then the mass of 
each star is determined by the integral, 



the volume occupied by each star is. 



[ p'dV, 
Jvi 



dV, 



and the spin angular momentum of each star is given by the expression. 



j: 



spm 



[ro^ sin^ + {wcosO - voif] Qp'dV . 



(30) 



(31) 



(32) 
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Having determined the volumes V* occupied by both rotationally flattened and tidally distorted 
stars, we define the mean radius of each star as, 

/oyi \ 1/3 

= ■ (^=> 

We furthermore define sub-domain T)^^ as the volume of the grid for which x > xli and 
^eff < ^e{i{xLi), and sub-domain as the volume of the grid for which x < xn and <I>cfr < 

^es{xLi)- Then the Roche-lobe volume surrounding each star is, 



RL 



/ dV, (34) 



and each star's Roche-lobe filling factor is obtained from the ratio, 



/rl = (35) 

■^RL 



3. Results 

3.1. Single White Dwarfs 

As mentioned earlier, we initially used a simplified version of our SCF code to construct a 
large number of single, nonrotating white dwarfs in order to compare our solutions with previous 
results (see Figure [T|) and to provide initial guesses for the density distributions inside both stars 
in each binary system. Table [2] details the properties of single, nonrotating white dwarfs that 
have central densities ranging from 10^'^ g cm~^ to lO^'' g cm~^ as determined from our model 
calculations; the 23 selected models are equally spaced in units of log Pmax- These spherical models 
were constructed on a uniform cylindrical mesh with resolution (128, 128, 128) in (ro, 6, z) using a 
convergence criterion 6 = 10~^. For each converged model, the first six columns of Table [2] list, 
respectively, the star's mass M in solar masses, radius R in units of 10® cm, central density Pmax 
in g cm~^, corresponding value of Xmax = (Pmax/ B)^^^ , moment of inertia, 

2 



I = J w^pdV , (36) 

in units of 10^'^ g cm^, and the radius of gyration, k = I /{ME?). As shown in the last column of 
Table [21 a typical virial error for these converged models was 10-"^ - 10-5. The values tabulated 
for the radius of gyration vary smoothly from k = 0.2036 for M = 0.0844Mq to = 0.1013 
for M = 1.4081Mq. This is consistent with our understanding that low- mass white dwarfs have 



structures similar to n = 3/2 polytropes for which k = 0.205 (jRucihskil Il988l ). while high-mas s 



white dwarfs display structures similar to n = 3 polytropes for which k = 0.0758 (|Rucihskilll988l ) 



Our values of k over this range of stellar masses are also consistent with the analytic function for 
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k{M) that iMarsh et al.l (|2004l ) fit through similar spherical model data. Knowledge of the radius 
of gyration of these spherical ZTWD models has assisted us in analyzing the tidally distorted 
structures that arise in our models of synchronously rotating white dwarfs in close binary systems 
(see further discussion, below). 

Using this same three-dimensional, cylindrical coordinate grid we constructed nonrotating 
models with central densities above lO^'^ g cm~^, that is, with masses above IAMq. We have not 
included these higher mass models in Table [2] or Figure [H however, because they did not converge 
to satisfactorily accurate structures. In particular, as the mass was steadily increased above IAMq, 
the models converged to structures with steadily increasing (rather than decreasing) values of k. 
By contrast, models constructed using a one-dimensional spherical code with much higher spatial 
resolution displayed values of k that decreased steadily to a value of 0.0755 at masses approaching 
Mch- If desired, the three-dimensional computational grid resolution could be increased to produce 
more accurate models of the white dwarf structure near the Chandrasekhar mass limit. 



3.2. White Dwarf Binary Sequences 

The slow inspiral evolution of a DWD binary can be mimicked by constructing a sequence of 
detached binaries having fixed Mtot and fixed q but varying separation, down to the separation 
at which the less massive star first makes contact with its Roche lobe. In an effort to illustrate 
the capabilities of our code, we have constructed three binary sequences having the same total 
mass — namely, Mtot = 1.5M0 — but three separate mass ratios. Specifically, sequence 'A' has 
g = 1, sequence 'B' has q = 2/3, and sequence 'C has q = 1/2. As detailed in Table EJ spherical 
models were constructed with the desired primary and secondary masses for these three sequences 
— specifically, M = 0.5Mq,0.6Mq,0.75Mq,0.9Mq and l.OM© — to provide good "guesses" for 
the initial binary star density distributions to start each SCF iteration. In addition to listing the 
values of M, R, Pmax, Xraa.^, and k for each of these converged spherical models, as was done for a 
wider range of spherical models in Table O Table [3] also lists values for the global energies W, U, 
and n in units of lO^'' ergs. 

Along each sequence, all the binary models were constructed using a uniform cylindrical grid 
with (128,256,65) zones in {w,6,z); by implementing reflection symmetry through the equatorial 
plane, only half as many zones were needed in the vertical direction as in the radial direction to 
achieve the same resolution in both. No additional symmetries were assumed in constructing the 
sequence, although, for the models shown here, symmetry through the x-z plane also could have 
been implemented for additional savings. The convergence criterion was set to 6 = 2.5 x 10~^; 
in most models, was the last variable to converge to this desired level. We note that, because 
the same number of grid zones was used for each model and each binary was scaled to fit entirely 
within the grid, the effective resolution of each star decreased as the binary separation a increased 
along each sequence. 
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Two tables have been produced for each DWD inspiral sequence in order to detail the properties 
of the models that lie along each sequence. For sequence 'A' {q = 1), for example, Table H] lists 
the values of six global binary system parameters (a, Q, Mtot, q, J tot, -E'tot) and the virial error 
obtained for thirty-five models (numbered Al through ^35) whose binary separations vary from 
2.0956 X 10^ cm at contact (model ^1) to 3.0911 x 10^ cm (model ^435). For this same group 
of models, Table [5] lists calculated values of five parameters (Mj, pJnaxi /rl' "^spin) 
individual stellar components (i = 1,2). Tables [6] and [7] provide the same detailed information for 
models along sequence 'B' {q = 2/3), and Tables [8] and [9] provide this information for models along 
sequence 'C [q = 1/2). 

The equatorial-plane density distributions displayed in Figures O HI and [5] illustrate the degree 
to which both white dwarf components are distorted by tides for various binary separations along 
each sequence. Labels in the upper-right-hand corner of each figure panel identify each binary 
system according to its corresponding position along each sequence, as itemized in Tables |4] - [9l 
Along sequence 'A' (Figure [3]), both components of the binary system are of equal size and display 
identical degrees of tidal distortion because the mass ratio q = 1. Along sequences 'B' and 'C 
(Figures U] and O respectively), however, the primary star (on the left in each figure panel) is 
noticeably smaller and less distorted than the secondary star. 

Figure [6] has been constructed from the data detailed in Tables [4] and [5] for binary sequence 
'A.' Specifically, the diamond symbols in the top two panels and in the bottom panel of this figure 
show, respectively, how the binary system's total angular momentum, Jtot, total energy, Etot-, and 
orbital angular velocity, $7, vary with orbital separation along this sequence; and the third panel 
from the top shows how the Roche-lobe filling factor, /fj^^, varies with orbital separation for both 
the primary star (diamonds) and the secondary star (asterisks). Figures [7| and [8] have been similarly 
constructed from the data detailed, respectively, in Tables [6] and [71 and in Tables [8] and [9l 



Following the lead of lNew fc Tohlind (119971 ). in constructing Figures [6] - [8] we have normalized 



our tabulated values of Jtot and E'tot to the quantities, 

Jnorm = (GM^. 75 -Ro.75 )'/' = 4.0735 X 10^° g cm^ s'S (37) 

Enorm = — ^ = 2.01 19 X 10^° erg , (38) 
-K0.75 

where i?o.75 = 7.4244 x 10^ cm is the radius of a spherical ZTWD whose mass is Mqj^ = 0.7522 Mq 
as tabulated in Table El Also, at each separation our tabulated values of have been normalized 
to the Keplerian orbital frequency. 

In all three figures, values of the orbital separation have been specified (bottom horizontal axis) in 
units of 10^ cm and (top horizontal axis) as normalized to the radius of a spherical ZTWD having 
the mass of the system's secondary star as tabulated in Table [3l that is, i2o.75 = 7.424 x 10^ cm, 
Rq qq = 8.671 X 10® cm and i^o.so = 9.638 x 10® cm. 
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New &: Tohlind (|l997l ) have previously constructed inspiral sequences for equal-mass DWD 
binary systems in which the structure of the individual component st ars is governed by the Chan- 
drasekhar ZTWD equation of state ([1]). The sequences published bv lNew &: Tohlind (|l997l ) cover 
a wide range of total masses. The one that most closely resembles our sequence 'A' (our only 
equal-mass sequence) has M^ot = l.G 3Mp,; the funct i onal b ehavior of £^tot(o) and Jtot(Q^) for this 
sequence is presented in Figure 16 of lNew &: Tohlind (Il997l ). Along this Mtot = 1.63M0 sequence, 



the two stars first make contact with their respective Roche lobes at a normalized separation of 
approximately 2.825 (see also Figure 5 of New & Tohline 1997). This is completely consistent with 
the behavior of our sequence 'A,' where contact occurs (model Al) when a/i?o.75 = 2.823. 



The DWD sequences constructed bv lNew &: Tohlind (jl997l ) all extend to separations smaller 
than the point of first contact, as their SCF technique allowed them to build over-contact (common 
envelope) equal-mass binaries. Their functions £'tot(a) and Jtot(«) display a quadratic behavior 
along the over-contact segment of each sequence, passing through a minimum at a binary separation 
smaller than the point of first contact. None of our three sequences show this behavior because we 
have not attempted to construct models past the initial point of contact. Indeed, it seems unlikely 
that equilibrium configurations exist at smaller separations except when the system mass ratio is 
precisely q = 1- 

For each of our DWD binary sequences, it is useful to compare the displayed functional behavior 
of ^tot(fl) from our numerical models against the behavior predicted by two simplified models. In 
the case of two point masses in circular orbit, the total angular momentum Jpm is given simply by 
the system's orbital angular momentum, that is, 



J, 



pm 



«^orb 



GMl,a 



(1 + 



(40) 



This function, normalized to Jnornn is displayed by the solid curve in the top panels of Figures [6] 
- El An even more realistic representation of the function Jtot(Q^) can be obtained by adding an 
approximate representation for the spin angular momentum, /jO, of both stars to the point-mass 
expression for Joi-b . If we assume that both stars retain a spherical structure while spinning at the 
Keplerian orbital frequency, Qk, the appropriate expression for the total "spinning sphere" system 
angular momentum is. 



J. 



Jorh + ih + h)^K = Jorb + [kiMiRt + fcsMsfl^ ) Qk 

2 / n \ 2n 



J, 



pm 



1 + 



1 + 



a 



+ 



R2 



(41) 



where, in addition to q, values of (the constants) Ri and ki appropriate for each binary sequence 
can be obtained from Table El Function (I4ip . normalized to Jnormi is displayed by the dot-dashed 
curve in the top panels of Figures[6l-El Analytic expression (j40p predicts that Jtot oc a^/^. Through 
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a correction factor, Eq. (j4ip displays a somewhat more complex behavior. Overall, our SCF model 
sequences match Eq. ()4ip particularly well. The largest deviation arises in all cases at the smallest 
separations; the slope of the SCF-generated Jtot(a) function flattens somewhat as the secondary 
star approaches contact with its Roche lobe, that is, as 1. 

The functional dependence of each system's total energy, Etot{a), can be understood in a 
similar fashion. Considering only the kinetic and gravitational potential energy of two point masses 
in circular orbit, we obtain. 



-Em- 



orb 



Korh + Worh = — -K'orb 



(42) 



where we have used the virial relation (2/Corb + W^orb) = 0. While this functional dependence 
explains the general £'tot(a) behavior exhibited in Figures [6] - [8] by our numerically constructed 
model sequences, expression (j42p is missing a nontrivial shift in the overall energy scale that is set 
by the binding energies of the two stars, namely. 



(43) 



Based on the properties of the spherical stellar models provided in Tabled the appropriate energy 
shift for sequences 'A,' 'B,' and 'C is, respectively, Ei, = —1.551 x 10^'^ ergs, —1.698 x 10^'^ ergs, 
and —1.963 x 10^*^ ergs. Adding E\y to i?orb provides what we will refer to as the "point mass" total 
system energy. 



E, 



pm 



+ ^b. 



(44) 



This analytic function, normalized to -Enormi is displayed as a solid curve in the plots of -Etot versus 
a shown in Figures [6] -[HI An improved approximation that we will refer to as the "spinning sphere" 
total system energy can be obtained by adding the rotational kinetic energy of both stars, assum- 
ing they remain spherically symmetric and spin uniformly with the Keplerian orbital frequency. 
Specifically, 



Ess — -Epm + "^ji^^i^K^ 

(1 + 



i=l 

Eh + Eorh { 1 



''It) M^) 



(45) 



This function, normalized to -Bnornn is displayed as a dot-dashed curve in the plots of E^tot versus 
a shown in Figures [6] - [HI Expression (j45p describes particularly well the variation of Stot with 
separation displayed by our numerically constructed binary sequences 'B' and 'C We note, however, 
that all three of our sequences show that the total system energy drops slightly below the behavior 
predicted by Eq. (j45p at the smallest separations. 
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The curve outlined by asterisks in the third panel from the top of Figures [6] - [8] shows that 
/rl^ steadily increases from a value ~ 0.2 to a value of 1.0 at the smallest separation along all three 
inspiral sequences, implying that the secondary star has made contact with its Roche lobe. For 
comparison, the curve outlined by diamonds in the same panel of these three figures shows how the 
Roche-lobe filling factor of the primary star varies along each sequence. The value of does not 
climb above 0.063 for sequence 'C or above 0.191 for sequence 'B,' reflecting the fact that in both 
cases the primary star is significantly more massive — and, hence, it has a significantly smaller 
radius — than the secondary star. For inspiral sequence 'A,' f^{a) displays an identical behavior 
to f^j^{cL) because the primary and secondary stars have equal masses. 

The bottom panel of Figures [6] - [8] displays the behavior of the normalized orbital frequency 
O/Ok as a function of binary separation derived from our three numerically constructed inspiral 
sequences. At the smallest separations, our models show that the orbital frequency is always 
~ 0.5% higher than predicted by the "point-mass" Keplerian frequency. Our equal-mass sequence 
exhibits the largest deviation at contact; specifically, for model 'Al,' we find Vt = 1.0085 Ok- As 
the separation is increased along each sequence, the figures show that O/^k approaches unity, 
as expected. However, at a sufficiently wide separation, each of our sequences displays a tiny 
discontinuous drop in the orbital frequency, followed by further decline that ultimately falls below 
the local Keplerian value. We suspect this odd behavior at wide separations arises from the discrete 
nature of our grid calculations coupled with progressively fewer grid zones falling inside both stars 
— resulting in progressively poorer numerical resolution — at wider separations. 



3.3. Conservative Mass- Transfer Sequences 

During a phase of stable mass transfer, a DWD binary system will evolve in such a way that the 
secondary star remains in marginal contact with its Roche lobe while it slowly transfers mass to the 
primary star. If the total mass of the system is conserved, then the evolution should proceed along 
a sequence of synchronously rotating configurations in which Mtot is constant, = 1, and q is 
steadily decreasing. Models Al, Bl and CI can be viewed as representing three such configurations 
along a sequence whose total system mass is Mtot = 1.5Mq. In evolving from an initially equal- 
mass, contact configuration (model ^41) to a semi-detached configuration with q = 2/3 (model -Bl), 
then on to a semi-detached configuration with q = 1/2 (model CI), the separation of such a system 
(measured in units of 10^ cm) will increase from og = 2.10, to ag = 2.67, then to ag = 3.18; and 
the system's orbital period (Porb = ^tt/Q) will increase from 42.3 s to 61.0 s, then to 79.5 s. 

It is clear, therefore, that our new SCF code can be used to construct model sequences that 
mimic the evolution of DWD systems undergoing slow, conservative mass-transfer. The models 
detailed in Tables [TOl (sequence 'D') and [TT] (sequence 'E') trace two such semi-detached sequences 
as the system mass ratio evolves from q = 1 to q < 0.5. For sequence 'D,' Mtot = 1.5Mq and for 
sequence 'E,' Mtot = l-OM©. In the top two panels of Figures [9l and [TOl data from Tables [TOl and 
[TT] have been plotted as diamond symbols to illustrate how a and Q vary with q while is held 
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to a value of unity (definition of a semi-detached binary) along these two fixed-mass sequences. 



Up to now, the community has relied upon some relatively simple analytic expressions to 
approximate the behavior of, for example, a{q) along conservative mass-transfer evolutionary tra- 
jectories. For example, by setting the radius of the secondary star as given by the Nauenberg 
mass-radius relation ()A14p equal to t he Roche- l obe r adius Rm, as defined in terms of a and q by 
the approximate relation provided by Eggleton ( 19831 ) . namely, 

0.4^/3 



Rrl 



0.6g2/3 + ln(l + gl/3) 



one obtains. 



where, 



^ ^ 0.0229(n2(z')~'/'(l 

Rq 



n 



4/3n1/2 



0.6g2/^ + ln(l + q 



1/3^ 



712 = 



(46) 



(47) 



(48) 



(1 + g) V Meh J 

The function a{q), defined by Eq. (|47|) for a given Mtot, has been plotted as a solid curve in the 
top panels of Figures and [TUl and the Keplerian orbital frequency associated with this separation 
(and relevant Mtot) has been plotted as a solid curve in the second panel of Figures [9l and [TOl For 
both sequence 'D' and sequence 'E,' the analytically derived curves are consistently offset by 3 - 
5% from our numerical model results. But overall, the analytically predicted functional behavior 
of a{q) and 0,{q) is in very good agreement with our results. This is reassuring as it provides a 
degree of validation for both our numerical code and the approximations that were adopted by 
earlier investigators when deriving the more easily manipulated analytic expressions. 

Finally, in the bottom two panels of Figures [9] and [TOl the diamond symbols display the 
variation of Jtot and E'tot with q along sequence 'D' and sequence 'E,' respectively. The solid 
curve drawn in the Jtot (9) panel of both figures shows the behavior predicted by our "spinning 
sphere" expression for the total system angular momentum (j4ip when used in conjunction with the 
o(g, Mtot) behavior prescribed by Eq. ([^7|) . Again, for a given Mtot there appears to be very good 
agreement between the functional behavior of Jtot(9) displayed by our numerical model results and 
the analytic expressions. There is also a systematic offset between the two. In either case it is clear 
that, unlike the behavior displayed by a{q) and ^{q), the system's total angular momentum does 
not vary monotonically with q along a conservative mass-transfer evolutionary trajectory. Note, in 
particular, that if the system mass ratio q is initially close to unity, Jtot increases as q decreases 
along the displayed trajectory. This result is unphysical. It signifies that slow evolution along 
a synchronously rotating, conservative mass-transfer trajectory can occur only if, at the onset of 
mass-transfer, q < gcrit, where the value of gcrit for a given Mtot is prescribed by the location of 
the maximum of the Jtotig) curve. For our model sequences 'D' and 'E,' we see that gcrit ^ 2/3, 
consistent with the mass-transfer stabi lity limit that has already received much attention in the 
literature ({Frank. Tohline k Evenlbood ). 
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4. Summary and Conclusions 



Based on the earlier work of tlie Hacliisu (1986a, b) and iHacIiisu. Eriguchi. &: Nomotd (|l986al ) 
we have developed a self-consistent-field technique that can be used to construct equilibrium models 
of synchronously rotating DWD binaries having a range of total masses, mass ratios, and binary 
separations. In addition to effects introduced by synchronous rotation, the distorted structure of 
both stars in each converged model is governed by the zero-temperature white dwarf equation of 
state ([T|) and a self-consistently determined, Newtonian gravitational field. In an effort to illustrate 
the technique's capabilities, we have constructed a set of models along five sequences: Three se- 
quences ('A', 'B', and 'C') mimic the last segment of the detached "inspiral" phase of DWD binary 
evolutions during which both -Mtot and q are held constant as a decreases; and two sequences ('D' 
and 'E') mimic a semi-detached "conservative mass transfer" phase of evolution during which Mtot 
is held fixed and the less massive star stays in marginal contact with its Roche lobe, but q steadily 
decreases while a steadily increases. 

Along each inspiral sequence, the functional dependence of Jtot and £^tot on the orbital sepa- 
ration can be well understood in terms of simple analytical expressions that describe two spinning 
spherical white dwarfs in circular orbit about one another. For a given total mass and separation, 
the calculated orbital frequencies along each inspiral sequence deviate measurably from associated 
Keplerian frequencies only in models for which the Roche-lobe filling factor of the less-massive star 
is > 60%. But, at least for the sequences examined here, the deviation from Keplerian frequencies 
is never more than 1% even at contact. 

Along both conservative mass-transfer sequences, we have documented how a, Jl, Jtot and E'tot 
vary with the system mass ratio as q decreases by roughly a factor of two, from q = 1.0 down 
to Q < 0.5. Along each sequence we have compared our numerically determined values of a at 
various values of q with the analytic a{q) function (H71) that is derived by setting the radius of the 



less massive star, as specified by the iNauenberd (|l972l ) mass-radius relation, equal to the Roche- 



lobe radius, as approximated bv lEggletonl (|l983l ). Qualitatively, our results show the same a{q) 



behavior that is predicted by this analytic expression. However, at a given q the value of a derived 
from our models is consistently ~ 8% larger than the value obtained from Eq. (j47p . The analytic 
expression could be brought into closer quantitative agreement with our numerical results if the 
leading coefficient in Eq. ()47p is increased by 8%, that is, if the expression's leading coefficient is 
changed from 0.0229 to 0.0247. This modification will, in turn, decrease the Keplerian frequency 
obtained from the analytic a{q) expression by ~ 9%, simultaneously bringing the analytically 
predicted orbital frequency into much closer agreement with our numerically determined values of 
Q. Along both of our conservative mass-transfer sequences, the plot of Jtotio) displays an extremum 
at a value of q < 2/3. The location of this extremum is almost certainly identifying the value of 
^crit that is relevant along both sequences. 

The development of this SCF technique was originally motivated by our desire to build models 
that would serve as good, "quiet" initial conditions for hydrodynamical simulations that are de- 
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signed to probe the onset and nonlinear development of mass-transfer instabilities in close, unequal- 
mass DWD binaries. The new computational tool that we have described in this paper achieves 
this objective. 



We acknowledge valuable interactions that we have had with B. Bourdin, J. Frank, D. Mar- 
cello, P. M. Motl, and S. Ou over the course of this project. We also thank an anonymous referee for 
pointing us to a key reference from Hachisu's collection of work during the mid-1980's. This work 
has been supported, in part, by grants AST-0708551 and DGE-0504507 from the U.S. National 
Science Foundation and, in part, by grant NNX07AG84G from NASA's ATP program. This re- 
search also has been made possible by grants of high-performance computing time on the TeraGrid 
(MCA98N043), at LSU, and across LONI (Louisiana Optical Network Initiative). 



A. White Dwarf Mass-Radius Relationship 
A.l. The Chandrasekhar Mass 



Chandrasekhaij (|l935l ) was the first to construct models of spherically symmetric stars using 
the equation of state defined by Eq. ([1]) and, in so doing, demonstrated that the maximum mass of 
an isolated, nonrotating white dwarf is M^h = l-44(//e/2)M0, where fie is the number of nucleons 
per electron and, hence, depends on the chemical composition of the WD. A conci se derivation of 
Mrh ( although, at the time, it was referred to as M3) is presented in Chapter XI of I Chandrasekhar 
(jl967l ). where we also find that the expressions for the two key coefficients in Eqs. ([l]) and ([2]) are. 



A 



4 5 
nm^c 

Sirrrir 



rrieC 



(Al) 
(A2) 



Numerical values for A and Bfi~^ are given here in Table lAll along with values of t he ph ysical 
constants c, h, rrie, and nip that we have used (column 2) and that IChandrasekhaii (|l967l ) used 
(column 3) to determine the values of A and Bfi~^. The derived analytic expression for the limiting 
mass is. 



2^y/^^ 

ttG) 



^ = 1.14205 X 10' 



,34 



where the coefficient. 



ma 



2.01824, 



(A3) 



(A4) 



represents a structural property of n = 3 polytrope s (7 = 4/3 gases) whose num erical value can 



be found in Chapter IV, Table 4 or IChandrasekharl (|l967l ). We note as well that [Chandrasekhar 
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(jl967l ) identified a characteristic radius, ii, for WDs given by the expression, 



= 7.71395 X 10^ cm . 

ttG / B 



(A5) 



A. 2. The "Nauenberg" Mass-Radius Relationship 



Nauenberel (|l972l ) derived an analytic approximation for the mass-radius relationship exhibited 



by isolated, spherical WDs that obey the ZTWD equation of state given in Eq. ([T]). Specifically, 
he offered an expression of the form. 
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(A6) 

(A7) 
(A8) 
(A9) 



is the atomic mass unit (see Table ET]) . fj, is the mean molecular weight of the gas, and C and i' 
are two adjustable parameters in Nauenberg's analytic approximation, both of which are expected 
to be of order unity. By assuming that the average particle mass denoted bv IChandrasekhaii (jl967l ) 
as {fierrip) is identical to the average particle mass specified by iNauenberd (jl972l ) as (fJ-m^) and, 
following Nauenberg's lead, by setting i/ = 1 ancj^. 



C = — = 0.324142 , 



in Eq. (jASp we see that, 



M, 



ch 



(AlO) 



(All) 



Hence, the denominator in (1A7P becomes the Chandrasekhar mass. Furthermore, expressions ()A9 
and ()A6P become, respectively. 



HeRo = m3(£i/ie) = 1.55686 X 10^ cm . 



(A12) 



Actually, iNauenberg j 19721 ) sets C = 0.323 
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and, 



R 



Rf 



[l-(M/Meh)"/^]^/=^ 

(M/A4h)i/3 



(A13) 



Finally, by adopting the values of Mq and Rq listed in Table lAll we obtain essentiall}o the 
identic al approxin i ate, a nalytic mass-radius relationship for ZTWDs presented in Eqs. (27) and 
(28) of lNauenber^ jmi ): 



R 

Rq 



0.0224 r [1 - (M/Mch)^/3]^/2 



(M/Meh)V3 



(A14) 



where. 



5.742 



(A15) 



A. 3. The "Eggleton" Mass-Radius Relationship 



Verbunt &: RappaportI (jl988l ) introduced the following approximate, analytic expression for 
the mass-radius relationship of a "completely degenerate . . . star composed of pure helium" (i.e.. 
He = 2), attributing its origin to Eggleton (private communication): 
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(A16) 



where Mp is a constant whose numerical value is O.OOO57M0 



This "Eggleton" mass-radius rela- 
tionship has been used widely by researchers when modeli ng the evolution o f semi-detached binary 
star systems in which the donor is a ZTWD. Since the iNauenbera (jl972l ) mass-radius relation- 
ship (|A14p is retrieved from Eq. ()A16P in the limit M/Mp ^ 1, it seems clear that Eggleton's 
contribution was th e insertion of the term in square brackets involving the ratio M/Adp which, as 



Marsh et al 



(120041') phrase it, "all ows for the change to a constant density configuration at low 



masses (jZapolskv fc Salpeterj|l969l )." In this paper we have only constructed binary star systems 
in which the internal structure of both stars is governed by the ZTWD equation of state ([T]). Hence 
it is appropriate for us to compare the properties of our modeled systems to behaviors predicted 
by the "Nauenberg," not the "Eggleton," mass-radius relationship. 



■^The numerical coefficients th at appear here in Eqs. (|A14|I and HA15|) differ slightly from the ones presented in 
Eqs. (27) and (28), respectively, of Nauenberg ( 19721 ) presumably because the values of the physical constants — such 
as Mq and Rq — that we have adopted in this paper (see Table IXTj) are slightly different from the values adopted 
by Nauenberg. 
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Table 1. Convergence of SCF Method: Binary Model B3 





Ne 






5 




VE 




64 


128 


33 


1.0 


X 10" 


-2 


4.5 X 10' 


-3 








1.0 X 10' 


-3 


2.6 X 10' 


-3 








1.4 X 10' 


-4 


2.2 X 10' 


-3 


128 


256 


65 


1.0 


X 10' 


-2 


4.0 X 10' 


-3 








1.0 


X 10' 


-3 


9.1 X 10' 


-4 








1.0 


X 10' 


-4 


5.7 X 10' 


-4 








3.5 


X 10' 


-5 


5.4 X 10' 


-4 
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Table 2. Sequence of single, nonrotating ZTWDs. 



M 
(Me) 


R 

(10^ cm) 


Pmax 

(g cm" 




^max 


/ 

(10^° g cm^) 




k 






VE 




n HQ /I A 


iy. (Dio 


3 1693 




iU 


U.zozo 




U. 


ZUOD 


1, 


.4 






-5 


n 1 1 1 o 
(J.iii<5 


17.9368 






10 


U.3060 


1 A A OO 

1.44ZZ 


0. 


ono 1 

,2031 


6, 


.1 






-5 


U.i4DU 


ib.ZDYZ 


1 onon 




10 


U.37U8 


l.ooUo 


0. 


ono A 

,z(Jz4 


1, 


.4 


V 




-5 


(J.i9U<5 


14.7421 






10 


U.44yz 


1.6517 


0, 


oni c; 

,2015 


5, 


O 

.5 


>^ 




-5 


(J. 2457 


13.3464 


3 1693 




10 


U.o44z 


1.7360 


0, 


onm 

,2001 


1, 


.4 


>^ 




-5 


(J.oio4 


12.0666 


0.6234 


X 


10 


0.6593 


1.7936 


0, 


,1983 


5, 


o 

.8 


X 


lU 


-5 


u.oyoo 


10.8906 


1.0000 


X 


iU 


0. (988 


1.8133 


0. 


,1958 


r 

0, 


.9 


X 


10" 


-5 


U.4o5y 


9.80/8 


1.7783 


X 


iU 


0.96 (8 


1. /8o2 


0. 


,1926 


-1 
i 


.6 


X 


10- 


-5 


0.5873 


8.8097 


3.1623 


X 


10*^ 


1.1725 


1.7051 


0, 


,1887 


6, 


.1 


X 


10- 


c: 
-0 


0.6942 


7.8896 


5.6234 


X 


10^ 


1.4205 


1.5754 


0, 


,1839 


5, 


.8 


X 


10- 


-5 


0.8018 


7.0419 


1.0000 


X 


10^ 


1.7209 


1.4058 


0, 


,1783 


7, 


.6 


X 


10- 


-5 


0.9058 


6.2628 


1.7783 


X 


10^ 


2.0850 


1.2122 


0, 


,1721 


6 


.6 


X 


10- 


-5 


1.0022 


5.5487 


3.1623 


X 


10^ 


2.5260 


1.0111 


0. 


,1653 


5, 


.1 


X 


10- 


-5 


1.0882 


4.8961 


5.6234 


X 


10^ 


3.0603 


0.8176 


0. 


,1581 


6, 


.3 


X 


10- 


-5 


1.1624 


4.3023 


1.0000 


X 


108 


3.7076 


0.6429 


0. 


,1507 


6 


.3 


X 


10- 


-5 


1.2246 


3.7643 


1.7783 


X 


108 


4.4919 


0.4929 


0. 


,1433 


7, 


.0 


X 


10- 


-5 


1.2753 


3.2793 


3.1623 


X 


108 


5.4421 


0.3698 


0. 


,1360 


7, 


.4 


X 


10- 


-5 


1.3155 


2.8443 


5.6234 


X 


108 


6.5932 


0.2723 


0. 


,1290 


9 


.0 


X 


10- 


-5 


1.3469 


2.4560 


1.0000 


X 


109 


7.9879 


0.1972 


0. 


,1225 


1, 


.1 


X 


10- 


-4 


1.3708 


2.1116 


1.7783 


X 


109 


9.6775 


0.1410 


0. 


,1164 


1, 


.2 


X 


10- 


-4 


1.3887 


1.8078 


3.1623 


X 


109 


11.7246 


0.0997 


0. 


,1108 


1, 


.3 


X 


10- 


-4 


1.4020 


1.5414 


5.6234 


X 


109 


14.2047 


0.0699 


0. 


,1058 


1, 


.4 


X 


10- 


-4 


1.4116 


1.3092 


1.0000 


X 


lOlo 


17.2094 


0.0486 


0. 


,1013 


1, 


.4 


X 


10- 


-4 



Table 3. Selected single, nonrotating ZTWDs. 



M 


R 


Pmax 




k 


W 


U 


n 


(Me) 


(10® cm) 


(g cm~'^) 






(10^° ergs) 




(10^° ergs) 


0.5019 


9.6383 


1.9536 X 10^ 


0.9979 


0.1923 


-0.6130 


0.3326 


0.2044 


0.6022 


8.6713 


3.4341 X 10^ 


1.2043 


0.1889 


-0.9932 


0.5526 


0.3311 


0.7522 


7.4244 


7.6648 X 10*^ 


1.5739 


0.1813 


-1.8512 


1.0757 


0.6171 


0.9028 


6.2847 


1.7483 X lO'^ 


2.0718 


0.1726 


-3.2433 


1.9860 


1.0812 


1.0025 


5.5456 


3.1703 X 10^ 


2.5265 


0.1655 


-4.6465 


2.9644 


1.5490 
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Table 4. DWD Inspiral Sequence 'A': Mtot = l.SM©; q = l 





n 


o 


-'"tot 




n 

y 


nnSO „ \ 


-E'tot 




VE 




Al 


2.0956 


14.8480 


1.5045 


1, 


.0000 


5.3879 


-1, 


.8624 


2.7 


X 


10- 


-4 


A2 


2.0970 


14.8317 


1.5043 


1, 


.0000 


5.3881 


-1, 


.8618 


2.8 


X 


10" 


-4 


A3 


2.1042 


14.7493 


1.5036 


1, 


.0000 


5.3882 


-1, 


.8589 


2.8 


X 


10" 


-4 


A4 


2.1099 


14.6847 


1.5030 


1, 


.0000 


5.3886 


-1, 


.8567 


2.9 


X 


10" 


-4 


A5 


2.1162 


14.6156 


1.5031 


1, 


.0000 


5.3915 


-1, 


.8566 


2.7 


X 


10" 


-4 


A6 


2.1239 


14.5339 


1.5031 


1, 


.0000 


5.3954 


-1, 


.8560 


2.7 


X 


10" 


-4 


A7 


2.1428 


14.3360 


1.5032 


1, 


.0000 


5.4059 


-1, 


.8549 


2.8 


X 


10" 


-4 


A8 


2.1544 


14.2154 


1.5030 


1, 


.0000 


5.4112 


-1, 


.8532 


2.9 


X 


10" 


-4 


A9 


2.1671 


14.0863 


1.5029 


1, 


.0000 


5.4180 


-1, 


.8519 


2.9 


X 


10" 


-4 


AlO 


2.1809 


13.9475 


1.5028 


1, 


.0000 


5.4254 


-1, 


.8508 


2.6 


X 


10" 


-4 


All 


2.1960 


13.7990 


1.5027 


1, 


.0000 


5.4337 


-1, 


.8490 


2.7 


X 


10" 


-4 


A12 


2.2292 


13.4849 


1.5030 


1, 


.0000 


5.4554 


-1, 


.8471 


2.6 


X 


10" 


-4 


A13 


2.2479 


13.3121 


1.5027 


1, 


.0000 


5.4665 


-1, 


.8448 


2.7 


X 


10" 


-4 


AU 


2.2669 


13.1429 


1.5032 


1, 


.0000 


5.4814 


-1, 


.8447 


2.7 


X 


10" 


-4 


A15 


2.2880 


12.9570 


1.5030 


1 


.0000 


5.4944 


-1, 


.8421 


2.7 


X 


10" 


-4 


^16 


2.3103 


12.7655 


1.5027 


1, 


.0000 


5.5082 


-1, 


.8392 


2.8 


X 


10" 


-4 


A17 


2.3572 


12.3804 


1.5029 


1 


.0000 


5.5423 


-1, 


.8359 


2.8 


X 


10" 


-4 


A18 


2.3828 


12.1772 


1.5027 


1, 


.0000 


5.5597 


-1, 


.8329 


2.8 


X 


10" 


-4 


A19 


2.4092 


11.9748 


1.5028 


1, 


.0000 


5.5792 


-1, 


.8307 


2.8 


X 


10" 


-4 


A20 


2.4362 


11.7751 


1.5032 


1, 


.0000 


5.6013 


-1, 


.8298 


2.8 


X 


10" 


-4 


A21 


2.4652 


11.5645 


1.5030 


1 


.0000 


5.6221 


-1, 


.8265 


2.9 


X 


10" 


-4 


A22 


2.5261 


11.1442 


1.5030 


1, 


.0000 


5.6684 


-1, 


.8215 


2.8 


X 


10" 


-4 


A23 


2.5584 


10.9314 


1.5030 


1, 


.0000 


5.6928 


-1, 


.8185 


2.9 


X 


10" 


-4 


A2A 


2.5920 


10.7176 


1.5028 


1, 


.0000 


5.7185 


-1, 


.8151 


3.1 


X 


10" 


-4 


A25 


2.6264 


10.5057 


1.5029 


1, 


.0000 


5.7457 


-1, 


.8127 


3.0 


X 


10" 


-4 


A26 


2.6624 


10.2904 


1.5028 


1, 


.0000 


5.7729 


-1, 


.8093 


3.0 


X 


10" 


-4 


A27 


2.6991 


10.0786 


1.5029 


1, 


.0000 


5.8016 


-1, 


.8068 


2.9 


X 


10" 


-4 


A28 


2.7376 


9.8598 


1.5028 


1, 


.0000 


5.8286 


-1, 


.8039 


1.5 


X 


10" 


-4 


A29 


2.8175 


9.4479 


1.5031 


1, 


.0000 


5.8981 


-1, 


.7984 


3.2 


X 


10" 


-4 


A30 


2.8597 


9.2377 


1.5031 


1, 


.0000 


5.9315 


-1, 


.7951 


3.2 


X 


10" 


-4 


A31 


2.8784 


9.1459 


1.5027 


1, 


.0000 


5.9440 


-1, 


.7925 


3.2 


X 


10" 


-4 


A32 


2.9478 


8.8242 


1.5032 


1, 


.0000 


6.0024 


-1, 


.7892 


3.3 


X 


10" 


-4 
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Table 4 — Continued 



Model 


a 




Mtot 




J tot 


Etot 


VE 




(10^ cm) 


(10-2 s-1) 


(Mo) 




(10^° cgs) 


(IqSO erg) 




^33 


2.9940 


8.6120 


1.5031 


1.0000 


6.0336 


-1.7860 


1.9 X 10"^ 


AM 


3.0420 


8.4102 


1.5031 


1.0000 


6.0738 


-1.7827 


2.1 X 10"^ 


^35 


3.0911 


8.2099 


1.5032 


1.0000 


6.1135 


-1.7797 


2.2 X lO"'' 



Table 5. Individual Stellar Components along DWD Inspiral Sequence 'A' 



Model Ml 


Ri 


i=i fi=i 


ji=l 
spin 


Af2 


R2 


i=2 fi=2 
Fmax J RL 


Ti=2 
spin 


(Mo) 


(10^ cm) 


(10^ cgs) 


(10^0 cgs) 


(Mo) 


(10^ cm) 


(10^ Cgs) 


(10^° cgs) 



Al 


0, 


.7522 


0. 


,7841 


6.745 


1, 


,0000 


0, 


.2562 


0, 


.7522 


0, 


.7841 


6.745 


1, 


,0000 


0, 


.2562 


A2 


0, 


.7522 


0. 


,7840 


6.745 


0, 


,9969 


0, 


.2559 


0, 


.7522 


0, 


.7840 


6.745 


0, 


.9969 


0, 


.2559 


A3 


0, 


.7518 


0. 


,7835 


6.745 


0, 


.9838 


0, 


.2540 


0, 


.7518 


0, 


.7835 


6.745 


0, 


.9838 


0, 


.2540 


A4 


0, 


.7515 


0. 


,7831 


6.745 


0, 


.9736 


0, 


.2526 


0, 


.7515 


0, 


.7831 


6.745 


0, 


.9736 


0, 


.2526 


A5 


0, 


.7516 


0. 


,7825 


6.758 


0, 


.9614 


0, 


.2509 


0, 


.7516 


0, 


.7825 


6.758 


0, 


.9615 


0, 


.2509 


A6 


0, 


.7516 


0. 


,7817 


6.771 


0, 


.9471 


0, 


.2489 


0, 


.7516 


0, 


.7817 


6.771 


0, 


.9471 


0, 


.2489 


A7 


0, 


.7516 


0. 


,7798 


6.805 


0, 


.9134 


0, 


.2441 


0, 


.7516 


0, 


.7798 


6.805 


0, 


.9134 


0, 


.2441 


A8 


0, 


.7515 


0. 


,7790 


6.820 


0, 


.8942 


0, 


.2414 


0, 


.7515 


0, 


.7790 


6.820 


0, 


.8942 


0, 


.2414 


A9 


0, 


.7514 


0. 


,7780 


6.837 


0, 


.8739 


0, 


.2384 


0, 


.7514 


0, 


.7780 


6.837 


0, 


.8739 


0, 


.2384 


AW 


0, 


.7514 


0. 


,7769 


6.855 


0, 


.8523 


0, 


.2353 


0, 


.7514 


0, 


.7769 


6.855 


0, 


.8523 


0, 


.2353 


All 


0, 


.7513 


0. 


,7759 


6.876 


0, 


.8303 


0, 


.2320 


0, 


.7513 


0, 


.7758 


6.876 


0, 


.8303 


0, 


.2320 


A12 


0, 


.7515 


0. 


,7735 


6.924 


0, 


.7838 


0, 


.2250 


0, 


.7515 


0, 


.7735 


6.924 


0, 


.7838 


0, 


.2250 


A13 


0, 


.7514 


0, 


,7724 


6.944 


0, 


.7602 


0, 


.2213 


0, 


.7514 


0, 


.7724 


6.944 


0, 


.7602 


0, 


.2213 


AU 


0, 


.7516 


0. 


,7711 


6.977 


0, 


.7367 


0, 


.2176 


0, 


.7516 


0, 


.7711 


6.977 


0, 


.7367 


0, 


.2176 


A15 


0, 


.7515 


0. 


,7701 


6.992 


0, 


.7131 


0, 


.2138 


0, 


.7515 


0, 


.7701 


6.992 


0, 


.7131 


0, 


.2138 


A16 


0, 


.7513 


0. 


,7692 


7.013 


0, 


.6903 


0, 


.2099 


0, 


.7513 


0, 


.7692 


7.013 


0, 


.6903 


0, 


.2099 


A17 


0, 


.7515 


0. 


,7668 


7.062 


0, 


.6440 


0, 


.2020 


0, 


.7515 


0, 


.7668 


7.062 


0, 


.6440 


0, 


.2020 


A18 


0, 


.7514 


0. 


,7658 


7.083 


0, 


.6215 


0, 


.1980 


0, 


.7514 


0, 


.7658 


7.083 


0, 


.6215 


0, 


.1980 


A19 


0, 


.7514 


0. 


,7650 


7.103 


0, 


.5998 


0, 


.1940 


0, 


.7514 


0, 


.7650 


7.103 


0, 


.5998 


0, 


.1940 


A20 


0, 


.7516 


0. 


,7637 


7.137 


0, 


.5779 


0, 


.1901 


0, 


.7516 


0, 


.7637 


7.137 


0, 


.5779 


0, 


.1901 


A21 


0, 


.7515 


0. 


,7629 


7.153 


0, 


.5568 


0, 


.1861 


0, 


.7515 


0, 


.7629 


7.153 


0, 


.5567 


0, 


.1861 


A22 


0, 


.7515 


0. 


,7611 


7.198 


0, 


.5160 


0, 


.1781 


0, 


.7515 


0, 


.7611 


7.198 


0, 


.5160 


0, 


.1781 


A23 


0, 


.7515 


0. 


,7603 


7.214 


0, 


.4964 


0, 


.1742 


0, 


.7515 


0, 


.7603 


7.214 


0, 


.4964 


0, 


.1742 


A24 


0, 


.7514 


0. 


,7596 


7.234 


0, 


.4774 


0, 


.1703 


0, 


.7514 


0, 


.7596 


7.234 


0, 


.4774 


0, 


.1703 



Table 5 — Continued 



Model Ml 


Ri 


t = l M = l 

Fmax ^RL 


p=l 

spin 


M2 


R2 


i=2 fi=2 


ji=2 
spin 


(Mo) 


(10^ cm) 


(10^ Cgs) 


(10^0 cgs) 


(Me) 


(10^ cm) 


(10^ Cgs) 


(10^° cgs) 



A25 


0, 


.7515 


0. 


,7585 


7.251 


0.4583 


0. 


,1664 


0, 


,7515 


0, 


.7585 


7.251 


0.4583 


0. 


,1664 


A26 


0, 


.7514 


0, 


,7579 


7.271 


0.4404 


0. 


,1625 


0, 


.7514 


0, 


.7579 


7.271 


0.4404 


0. 


,1625 


A27 


0, 


.7514 


0, 


,7571 


7.290 


0.4231 


0. 


,1587 


0, 


.7514 


0, 


.7571 


7.290 


0.4231 


0. 


,1587 


A28 


0, 


.7514 


0. 


,7562 


7.309 


0.4057 


0. 


,1548 


0, 


.7514 


0, 


.7562 


7.309 


0.4057 


0. 


,1548 


A29 


0, 


.7516 


0. 


,7548 


7.349 


0.3733 


0. 


,1476 


0, 


.7516 


0, 


.7548 


7.349 


0.3733 


0. 


,1476 


A30 


0, 


.7515 


0. 


,7541 


7.363 


0.3575 


0. 


,1440 


0, 


.7515 


0, 


.7541 


7.363 


0.3576 


0. 


,1440 


A31 


0, 


.7513 


0, 


,7540 


7.363 


0.3478 


0. 


,1424 


0, 


.7513 


0, 


.7541 


7.363 


0.3518 


0. 


,1424 


A32 


0, 


.7516 


0. 


,7529 


7.397 


0.3282 


0. 


,1369 


0, 


.7516 


0, 


.7529 


7.397 


0.3282 


0. 


,1369 


A33 


0, 


.7515 


0. 


,7522 


7.412 


0.3142 


0. 


,1333 


0, 


.7515 


0, 


.7522 


7.412 


0.3142 


0. 


,1333 


A34 


0, 


.7516 


0. 


,7517 


7.425 


0.3005 


0. 


,1299 


0, 


.7516 


0, 


.7517 


7.425 


0.3005 


0. 


,1299 


A35 


0, 


.7516 


0. 


,7511 


7.438 


0.2874 


0. 


1265 


0, 


.7516 


0, 


.7510 


7.440 


0.2873 


0. 


,1265 
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Table 6. DWD Inspiral Sequence 'B': Mtot = l.SM©; q = 2/3 



A/Tnrlpl 

iVlUU-Ci 


a 

lU Clll j 


o 

[W S ) 


A/L ^ 


Q 


'J tot 
[lU CgS) 


-E'tot 




VE 




Bl 


2.6679 


10.2944 


1.5042 


0.6671 


5.5888 


-1, 


.9460 


6.0 


X 


10" 


-4 


B2 


2.6743 


10.2576 


1.5043 


0.6667 


5.5931 


-1, 


.9460 


6.1 


X 


10- 


-4 


B3 


2.6819 


10.2106 


1.5037 


0.6667 


5.5951 


-1, 


.9436 


6.0 


X 


10- 


-4 


B4 


2.6923 


10.1491 


1.5033 


0.6664 


5.6001 


-1, 


.9419 


6.1 


X 


10- 


-4 


B5 


2.7039 


10.0831 


1.5036 


0.6663 


5.6090 


-1, 


.9420 


6.0 


X 


10- 


-4 


B6 


2.7177 


10.0040 


1.5034 


0.6665 


5.6174 


-1, 


.9401 


6.1 


X 


10- 


-4 


B7 


2.7487 


9.8330 


1.5035 


0.6665 


5.6403 


-1, 


.9383 


6.2 


X 


10- 


-4 


B8 


2.7687 


9.7249 


1.5035 


0.6664 


5.6542 


-1, 


.9372 


6.2 


X 


10- 


-4 


B9 


2.7902 


9.6104 


1.5033 


0.6665 


5.6684 


-1, 


.9348 


6.3 


X 


10- 


-4 


BIO 


2.8138 


9.4879 


1.5033 


0.6664 


5.6851 


-1, 


.9333 


6.3 


X 


10- 


-4 


Bll 


2.8388 


9.3609 


1.5033 


0.6663 


5.7034 


-1, 


.9318 


6.4 


X 


10- 


-4 


B12 


2.8654 


9.2290 


1.5032 


0.6665 


5.7224 


-1, 


.9295 


6.4 


X 


10- 


-4 


B13 


2.9196 


8.9712 


1.5034 


0.6665 


5.7638 


-1, 


.9267 


6.5 


X 


10- 


-4 


BU 


2.9514 


8.8246 


1.5033 


0.6665 


5.7864 


-1, 


.9241 


6.8 


X 


10- 


-4 


B15 


2.9850 


8.6746 


1.5031 


0.6665 


5.8108 


-1, 


.9215 


6.8 


X 


10- 


-4 


B16 


3.0199 


8.5240 


1.5034 


0.6664 


5.8380 


-1, 


.9201 


6.9 


X 


10- 


-4 


B17 


3.0567 


8.3680 


1.5031 


0.6665 


5.8638 


-1, 


.9170 


7.0 


X 


10- 


-4 


B18 


3.0949 


8.2123 


1.5032 


0.6665 


5.8928 


-1, 


.9147 


7.0 


X 


10- 


-4 


B19 


3.1350 


8.0548 


1.5034 


0.6663 


5.9238 


-1, 


.9132 


7.2 


X 


10- 


-4 


B20 


3.1766 


7.8955 


1.5031 


0.6665 


5.9543 


-1, 


.9097 


7.4 


X 


10- 


-4 


B21 


3.2199 


7.7355 


1.5032 


0.6666 


5.9870 


-1, 


.9074 


7.5 


X 


10- 


-4 


B22 


3 2653 


7.5743 


1.5034 


0.6664 


6.0218 


-1, 


.9056 


7.8 


X 


10- 


-4 


B23 


3.3554 


7.2699 


1.5033 


0.6667 


6.0905 


-1, 


.9002 


8.1 


X 


10- 


-4 


B24: 


3.4062 


7.1070 


1.5034 


0.6666 


6.1285 


-1, 


.8978 


8.1 


X 


10- 


-4 


B25 


3.4590 


6.9440 


1.5034 


0.6665 


6.1681 


-1, 


.8951 


8.5 


X 


10- 


-4 


B2% 


3.5133 


6.7821 


1.5032 


0.6666 


6.2078 


-1, 


.8917 


8.6 


X 


10- 


-4 


B27 


3.5700 


6.6211 


1.5033 


0.6665 


6.2506 


-1, 


.8893 


8.8 


X 


10- 


-4 


B28 


3.6285 


6.4611 


1.5034 


0.6665 


6.2948 


-1, 


.8867 


9.2 


X 


10- 


-4 


B29 


3.6892 


6.3010 


1.5032 


0.6666 


6.3387 


-1, 


.8832 


9.3 


X 


10- 


-4 


B30 


3.7520 


6.1432 


1.5033 


0.6666 


6.3860 


-1, 


.8809 


9.6 


X 


10- 


-4 


B31 


3.8171 


5.9860 


1.5034 


0.6665 


6.4340 


-1, 


.8783 


9.9 


X 


10- 


-4 


B32 


3.8846 


5.8298 


1.5032 


0.6667 


6.4828 


-1, 


.8747 


1.0 


X 


10- 


-3 
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Table 6 — Continued 



Model 


a 


n 


A^tot 


Q 




-E'tot 






VE 






(10^ cm) 


/ . „ O 1 V 

(10-2 s-^) 






(10^° cgs) 


(10^0 


erg) 












B33 


3.9544 


5.6759 


1.5034 


0.6666 


6.5344 


-1, 


.8724 


1, 


.0 


X 


10" 


-3 


B34 


4.0270 


5.5225 


1.5034 


0.6665 


6.5874 


-1, 


.8698 


1, 


.1 


X 


10' 


-3 


535 


4.1030 


5.3681 


1.5027 


0.6666 


6.6380 


-1, 


.8645 


1, 


.1 


X 


10' 


-3 


B36 


4.1801 


5.2199 


1.5031 


0.6666 


6.6951 


-1, 


.8628 


1, 


.1 


X 


10' 


-3 


B37 


4.2601 


5.0741 


1.5034 


0.6666 


6.7548 


-1, 


.8611 


1, 


.2 


X 


10' 


-3 


B38 


4.2667 


5.0582 


1.5030 


0.6668 


6.7529 


-1, 


.8596 


1, 


.1 


X 


10' 


~3 


B39 


4.3502 


4.9129 


1.5031 


0.6667 


6.8124 


-1, 


.8570 


1, 


.2 


X 


10' 


-3 


540 


4.4375 


4.7683 


1.5032 


0.6664 


6.8739 


-1, 


.8548 


1, 


.2 


X 


10' 


-3 


541 


4.5270 


4.6269 


1.5030 


0.6666 


6.9357 


-1, 


.8513 


1, 


.3 


X 


10' 


-3 


542 


4.6192 


4.4889 


1.5032 


0.6668 


7.0011 


-1, 


.8487 


1, 


.3 


X 


10' 


-3 


543 


4.7160 


4.3508 


1.5034 


0.6666 


7.0680 


-1, 


.8467 


1, 


.4 


X 


10" 


-3 


544 


4.8164 


4.2143 


1.5031 


0.6667 


7.1344 


-1, 


.8429 


1, 


.4 


X 


10" 


-3 


545 


4.9203 


4.0808 


1.5029 


0.6668 


7.2039 


-1, 


.8395 


1, 


.5 


X 


10" 


-3 


546 


5.0281 


3.9494 


1.5029 


0.6668 


7.2752 


-1, 


.8366 


1, 


.5 


X 


10" 


-3 


547 


5.1405 


3.8227 


1.5033 


0.6665 


7.3562 


-1, 


.8352 


1, 


.7 


X 


10" 


-3 


548 


5.1467 


3.8136 


1.5032 


0.6669 


7.3567 


-1, 


.8347 


1, 


.6 


X 


10" 


-3 


549 


5.2630 


3.6870 


1.5033 


0.6670 


7.4330 


-1, 


.8320 


1, 


.7 


X 


10" 


-3 



Table 7. Individual Stellar Components along DWD Inspiral Sequence 'B' 



Model Ml R^ p-^ 4=i M2 R2 f^=l f^L «^sV1n 





(Mo) 


(10^ cm) 


(10^ cgs) 






(10^0 


cgs) 


(Mo) 


(10^ cm) 


(109 . 


cgs) 






(1050 


cgs) 


Bl 


0.9023 


0.6345 


17.049 


0. 


1913 


0. 


,1289 


0.6019 


0.9117 


3, 


.076 


1, 


.0001 


0. 


,2004 


B2 


0.9026 


0.6344 


17.049 


0. 


1898 


0. 


,1283 


0.6017 


0.9113 


3, 


.076 


0, 


.9908 


0. 


,1995 


B3 


0.9022 


0.6345 


17.049 


0. 


1881 


0. 


,1278 


0.6015 


0.9108 


3, 


.076 


0, 


.9805 


0. 


,1982 


BA 


0.9022 


0.6343 


17.049 


0. 


1856 


0. 


,1270 


0.6012 


0.9101 


3, 


.076 


0, 


,9667 


0. 


,1966 


B5 


0.9023 


0.6343 


17.047 


0. 


1832 


0. 


,1261 


0.6013 


0.9090 


3, 


.083 


0, 


,9504 


0. 


,1948 


B6 


0.9021 


0.6343 


17.059 


0. 


1803 


0. 


,1250 


0.6013 


0.9078 


3, 


,089 


0, 


.9314 


0. 


,1927 


B7 


0.9022 


0.6340 


17.082 


0. 


1737 


0. 


,1228 


0.6013 


0.9056 


3, 


,108 


0, 


,8916 


0. 


,1882 


B8 


0.9022 


0.6336 


17.085 


0. 


1694 


0. 


,1213 


0.6012 


0.9043 


3, 


.115 


0, 


.8677 


0. 


,1854 


B9 


0.9020 


0.6336 


17.100 


0. 


1654 


0. 


,1198 


0.6012 


0.9029 


3, 


,123 


0, 


.8430 


0. 


,1825 


BW 


0.9021 


0.6334 


17.112 


0. 


1610 


0. 


,1182 


0.6011 


0.9017 


3, 


,134 


0, 


.8181 


0. 


,1795 


Bll 


0.9022 


0.6334 


17.121 


0. 


1565 


0. 


,1165 


0.6011 


0.9003 


3, 


,144 


0, 


.7917 


0. 


,1764 


B12 


0.9020 


0.6332 


17.139 


0. 


1520 


0. 


,1148 


0.6012 


0.8987 


3, 


,153 


0, 


.7650 


0. 


,1732 


B13 


0.9022 


0.6329 


17.162 


0. 


1432 


0. 


,1114 


0.6013 


0.8963 


3, 


,175 


0, 


.7172 


0. 


,1670 


BU 


0.9021 


0.6326 


17.177 


0. 


1383 


0. 


,1095 


0.6012 


0.8950 


3, 


.184 


0, 


,6907 


0. 


,1636 


B15 


0.9020 


0.6324 


17.195 


0. 


1335 


0. 


,1076 


0.6011 


0.8935 


3, 


.194 


0, 


,6644 


0. 


,1602 


B16 


0.9022 


0.6323 


17.203 


0. 


1288 


0. 


,1056 


0.6012 


0.8925 


3, 


.203 


0, 


,6403 


0. 


,1568 


B17 


0.9020 


0.6323 


17.222 


0. 


1241 


0. 


1036 


0.6011 


0.8914 


3, 


.213 


0, 


,6158 


0. 


,1534 


B18 


0.9020 


0.6318 


17.241 


0. 


1192 


0. 


,1016 


0.6012 


0.8902 


3, 


.223 


0, 


.5916 


0. 


,1500 


B19 


0.9022 


0.6316 


17.246 


0. 


1145 


0. 


,0995 


0.6012 


0.8887 


3, 


.233 


0, 


.5675 


0. 


,1465 


B20 


0.9019 


0.6315 


17.268 


0. 


1099 


0. 


,0975 


0.6012 


0.8879 


3, 


.242 


0, 


.5448 


0. 


,1431 


B21 


0.9020 


0.6311 


17.287 


0. 


1053 


0. 


,0955 


0.6012 


0.8866 


3, 


.251 


0, 


.5219 


0. 


,1397 


B22 


0.9022 


0.6311 


17.294 


0. 


1009 


0. 


,0933 


0.6012 


0.8856 


3, 


,261 


0, 


,5003 


0. 


,1364 


B23 


0.9020 


0.6307 


17.340 


0. 


,0928 


0. 


,0895 


0.6013 


0.8837 


3, 


.280 


0, 


,4605 


0. 


,1301 


B24 


0.9021 


0.6304 


17.353 


0. 


,0885 


0. 


,0874 


0.6013 


0.8829 


3, 


.285 


0, 


,4408 


0. 


,1268 



Table 7 — Continued 



Model Ml 


^1 rmax 


fi=l 
JRL 


Ti=l 

spm 


M2 


^2 rmax 


fi=2 
JRL 


Ti=2 
spin 


(Mo) 


(10^ cm) (10^ cgs) 




(10^° cgs) 


(Mo) 


(10^ cm) (10^ cgs) 




(10^° cgs) 



B25 


0, 


.9021 


0, 


.6302 


17.362 


0, 


,0844 


0.0853 


0, 


,6012 


0.8823 


3, 


.294 


0, 


,4219 


0, 


,1235 


B26 


0, 


.9019 


0, 


.6303 


17.384 


0, 


.0806 


0.0833 


0, 


.6013 


0.8814 


3, 


.301 


0, 


,4029 


0, 


,1203 


B27 


0, 


.9021 


0, 


.6301 


17.398 


0, 


.0767 


0.0812 


0, 


,6012 


0.8804 


3, 


.309 


0, 


,3846 


0, 


,1171 


B28 


0, 


.9021 


0. 


.6297 


17.409 


0, 


.0729 


0.0792 


0, 


,6013 


0.8791 


3, 


.316 


0, 


,3661 


0, 


,1140 


B29 


0, 


.9019 


0. 


.6298 


17.432 


0, 


.0694 


0.0772 


0, 


,6013 


0.8790 


3, 


.323 


0, 


,3500 


0, 


,1109 


B30 


0, 


.9021 


0. 


.6292 


17.445 


0, 


.0657 


0.0752 


0, 


.6013 


0.8780 


3, 


.330 


0, 


,3334 


0, 


,1078 


B31 


0, 


.9021 


0. 


.6293 


17.457 


0, 


.0624 


0.0732 


0, 


.6013 


0.8772 


3, 


.336 


0, 


,3174 


0, 


,1048 


B32 


0, 


.9019 


0, 


.6288 


17.480 


0, 


.0591 


0.0712 


0, 


.6013 


0.8768 


3, 


.343 


0, 


,3024 


0, 


,1019 


i?33 


0, 


.9020 


0, 


.6289 


17.495 


0, 


.0560 


0.0693 


0, 


.6013 


0.8762 


3, 


.349 


0, 


,2879 


0, 


,0990 


534 


0, 


.9021 


0, 


.6291 


17.504 


0, 


.0531 


0.0674 


0, 


.6013 


0.8753 


3, 


.355 


0, 


,2734 


0, 


,0961 


B35 


0, 


.9017 


0, 


.6287 


17.504 


0, 


.0501 


0.0655 


0, 


.6011 


0.8749 


3, 


.355 


0, 


,2597 


0, 


,0932 


i?36 


0, 


.9019 


0, 


.6285 


17.531 


0, 


.0473 


0.0636 


0, 


.6012 


0.8746 


3, 


.365 


0, 


,2469 


0, 


,0905 


B37 


0, 


.9021 


0, 


.6283 


17.542 


0, 


.0447 


0.0617 


0, 


.6013 


0.8739 


3, 


.371 


0, 


,2344 


0, 


,0878 


B38 


0, 


.9017 


0, 


.6284 


17.542 


0, 


.0445 


0.0616 


0, 


.6013 


0.8739 


3, 


.371 


0, 


,2334 


0, 


,0875 


B39 


0, 


.9018 


0, 


.6280 


17.572 


0, 


.0419 


0.0598 


0, 


.6013 


0.8732 


3, 


.378 


0, 


,2212 


0, 


,0848 


540 


0, 


.9021 


0, 


.6276 


17.572 


0, 


.0394 


0.0579 


0, 


.6011 


0.8728 


3, 


.378 


0, 


,2096 


0, 


,0822 


541 


0, 


.9018 


0, 


.6273 


17.601 


0, 


.0371 


0.0562 


0, 


.6012 


0.8723 


3, 


.387 


0, 


,1984 


0, 


,0796 


542 


0, 


.9018 


0, 


.6274 


17.631 


0, 


.0350 


0.0545 


0, 


.6014 


0.8712 


3, 


.392 


0, 


,1875 


0, 


,0771 


543 


0, 


.9021 


0, 


.6267 


17.645 


0, 


.0328 


0.0527 


0, 


.6013 


0.8710 


3, 


.398 


0, 


,1774 


0, 


,0746 


544 


0, 


.9019 


0, 


.6272 


17.645 


0, 


.0309 


0.0510 


0, 


.6012 


0.8709 


3, 


.398 


0, 


,1677 


0, 


,0722 


545 


0, 


.9017 


0, 


.6272 


17.682 


0, 


.0290 


0.0494 


0, 


.6012 


0.8704 


3, 


.407 


0, 


,1583 


0, 


,0698 


546 


0, 


.9017 


0, 


.6265 


17.682 


0, 


.0271 


0.0478 


0, 


.6012 


0.8699 


3, 


.407 


0, 


,1491 


0, 


,0674 


547 


0, 


.9020 


0, 


.6266 


17.697 


0, 


.0254 


0.0462 


0, 


.6012 


0.8696 


3, 


.416 


0, 


,1405 


0, 


,0652 


548 


0, 


.9018 


0, 


.6266 


17.726 


0, 


.0254 


0.0461 


0, 


.6014 


0.8695 


3, 


.419 


0, 


,1400 


0, 


,0650 



Table 7— Continued 



Model 


Ml 
(Mo) 


Ri 
(10^ cm) 


r max 

(10^ cgs) 


fi=l 
JRL 


ji=i 

spm 

(10^° cgs) 


M2 

(Mo) 


R2 

(10^ cm) 


p'=2 

r max 

(10^ cgs) 


fi=2 
JRL 


ji=2 
spin 

(10^° cgs) 


BA9 


0.9018 


0.6260 


17.758 


0.0237 


0.0445 


0.6015 


0.8691 


3.421 


0.1318 


0.0628 
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Table 8. DWD Inspiral Sequence 'C: Mtot = l-SM©; q = l/2 



A/Tnrlpl 

iVlUU-Ci 


a 


o 

s ) 


-'Wtot 


Q 


'J tot 
[lU CgS) 


-E'tot 






VE 




CI 


3.1807 


7.9054 


1.5043 


0.5001 


5.5568 


-2. 


,1635 


1, 


.2 


X 


10" 


-3 


C2 


3.1814 


7.9025 


1.5044 


0.5001 


5.5574 


-2. 


,1638 


1, 


.2 


X 


10- 


-3 


C3 


3.1865 


7.8822 


1.5039 


0.5002 


5.5582 


-2. 


,1615 


1, 


.2 


X 


10- 


-3 


C4 


3.2070 


7.8041 


1.5033 


0.4997 


5.5671 


-2. 


,1594 


1, 


.2 


X 


10- 


-3 


C5 


3.2208 


7.7539 


1.5038 


0.4997 


5.5780 


-2. 


,1605 


1, 


.2 


X 


10- 


-3 


C6 


3.2376 


7.6922 


1.5036 


0.4998 


5.5883 


-2. 


,1584 


1, 


.2 


X 


10- 


-3 


C7 


3.2573 


7.6213 


1.5034 


0.4999 


5.6009 


-2. 


,1566 


1, 


.2 


X 


10- 


-3 


C8 


3.2800 


7.5404 


1.5035 


0.4998 


5.6150 


-2. 


,1561 


1, 


.2 


X 


10- 


-3 


C9 


3.3051 


7.4547 


1.5038 


0.4996 


5.6329 


-2. 


,1566 


1, 


.2 


X 


10- 


-3 


CIO 


3.3321 


7.3618 


1.5033 


0.4999 


5.6489 


-2. 


,1529 


1, 


.2 


X 


10- 


-3 


Cll 


3.3617 


7.2631 


1.5032 


0.4999 


5.6680 


-2. 


,1511 


1, 


.2 


X 


10- 


-3 


C12 


3.4235 


7.0659 


1.5036 


0.4999 


5.7108 


-2. 


,1497 


1, 


.3 


X 


10- 


-3 


CIS 


3.4601 


6.9530 


1.5036 


0.4998 


5.7349 


-2. 


,1483 


1, 


.3 


X 


10- 


-3 


C14 


3.4985 


6.8369 


1.5032 


0.4999 


5.7590 


-2. 


,1448 


1, 


.3 


X 


10- 


-3 


C15 


3.5391 


6.7190 


1.5033 


0.5000 


5.7867 


-2. 


,1432 


1, 


.3 


X 


10- 


-3 


C16 


3.5823 


6.5967 


1.5035 


0.4998 


5.8157 


-2. 


,1426 


1, 


.4 


X 


10- 


-3 


C17 


3.6274 


6.4733 


1.5035 


0.4998 


5.8463 


-2. 


,1405 


1, 


.4 


X 


10- 


-3 


CIS 


3.6747 


6.3467 


1.5032 


0.5000 


5.8765 


-2. 


,1370 


1, 


.4 


X 


10- 


-3 


C19 


3.7245 


6.2193 


1.5032 


0.4999 


5.9102 


-2. 


,1353 


1, 


.4 


X 


10- 


-3 


C20 


3.7764 


6.0914 


1.5036 


0.4998 


5.9468 


-2. 


,1349 


1, 


.5 


X 


10- 


-3 


C21 


3.8311 


5.9605 


1.5034 


0.4998 


5.9826 


-2. 


,1320 


1, 


.5 


X 


10- 


-3 


C22 


3.8878 


5.8291 


1.5032 


0.5001 


6.0197 


-2. 


,1286 


1, 


.5 


X 


10- 


-3 


C23 


3.9475 


5.6967 


1.5032 


0.4999 


6.0592 


-2. 


,1266 


1, 


.6 


X 


10- 


-3 


C24 


4.0093 


5.5658 


1.5036 


0.4998 


6.1025 


-2. 


,1261 


1, 


.6 


X 


10- 


-3 


C25 


4.0742 


5.4324 


1.5034 


0.4999 


6.1443 


-2. 


,1232 


1, 


.7 


X 


10- 


-3 


C26 


4.1412 


5.2995 


1.5031 


0.5001 


6.1874 


-2. 


,1194 


1, 


.7 


X 


10- 


-3 


C27 


4.2115 


5.1674 


1.5031 


0.5000 


6.2340 


-2. 


,1173 


1, 


.8 


X 


10- 


-3 


C28 


4.2844 


5.0359 


1.5035 


0.4999 


6.2833 


-2. 


,1166 


1, 


.8 


X 


10- 


-3 


C29 


4.3605 


4.9039 


1.5035 


0.4999 


6.3327 


-2. 


,1146 


1, 


.9 


X 


10- 


-3 


C30 


4.4396 


4.7723 


1.5031 


0.5001 


6.3821 


-2. 


,1100 


1, 


.9 


X 


10- 


-3 


C31 


4.5218 


4.6425 


1.5031 


0.5001 


6.4355 


-2. 


,1076 


2, 


.0 


X 


10- 


-3 


C32 


4.6077 


4.5130 


1.5033 


0.5000 


6.4912 


-2. 


,1066 


2, 


.0 


X 


10- 


-3 
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Table 8 — Continued 



Model 


a 


n 


A^tot 


q 


Jtot 


-E'tot 


VE 






(10^ cm) 




(Me) 




(10^° cgs) 


(IqSO gj.g) 






C33 


4.6142 


4.5018 


1.5040 


0.4996 


6.4939 


-2.1098 


2.0 X 10' 


-3 


C34 


4.7030 


4.3735 


1.5033 


0.5001 


6.5480 


-2.1038 


2.1 X 10" 


-3 


C35 


4.7960 


4.2462 


1.5030 


0.5001 


6.6056 


-2.1003 


2.2 X 10" 


-3 


C36 


4.8938 


4.1186 


1.5027 


0.5000 


6.6642 


-2.0972 


2.2 X 10" 


-3 



Table 9. Individual Stellar Components along DWD Inspiral Sequence 'C 



Model Ml 


Ri 


i=i fi=i 


ji=l 
spin 


Af2 


R2 


i=2 fi=2 
Fmax J RL 


Ti=2 
spin 


(Mo) 


(10^ cm) 


(10^ Cgs) 


(10^0 cgs) 


(Mo) 


(10^ cm) 


(10^ Cgs) 


(10^° cgs) 



CI 


1 


.0028 


0. 


,5542 


31, 


.824 


0.0627 


0.0800 


0, 


,5015 


1 


.0110 


1 


.767 


1, 


,0000 


0. 


,1613 


C2 


1 


.0029 


0. 


,5541 


31 


.824 


0.0626 


0.0799 


0, 


,5015 


1 


.0108 


1 


.767 


0, 


.9987 


0. 


,1613 


C3 


1 


.0025 


0. 


,5544 


31, 


.824 


0.0624 


0.0798 


0, 


.5014 


1 


.0104 


1 


.767 


0, 


,9927 


0. 


,1607 


C4 


1 


.0025 


0. 


,5546 


31 


.824 


0.0612 


0.0790 


0, 


.5009 


1 


.0093 


1 


.767 


0, 


.9703 


0. 


,1586 


C5 


1, 


.0027 


0. 


,5541 


31, 


.806 


0.0602 


0.0784 


0, 


,5010 


1, 


.0081 


1 


.772 


0, 


.9536 


0. 


,1571 


C6 


1 


.0025 


0. 


,5543 


31, 


.812 


0.0593 


0.0778 


0, 


.5011 


1 


.0068 


1 


.776 


0, 


.9334 


0. 


,1553 


C7 


1 


.0023 


0, 


,5547 


31 


.834 


0.0583 


0.0771 


0, 


.5011 


1 


.0051 


1 


.780 


0, 


.9109 


0. 


,1534 


C8 


1 


.0024 


0. 


,5541 


31, 


.845 


0.0568 


0.0762 


0, 


,5010 


1 


.0035 


1 


.786 


0, 


.8873 


0. 


,1512 


C9 


1 


.0028 


0. 


,5540 


31, 


.823 


0.0554 


0.0753 


0, 


,5010 


1 


.0021 


1 


.791 


0, 


.8634 


0. 


,1489 


CIO 


1 


.0023 


0. 


,5544 


31 


.847 


0.0542 


0.0744 


0, 


.5010 


1, 


.0007 


1 


.796 


0, 


.8375 


0. 


,1464 


Cll 


1 


.0022 


0. 


,5541 


31 


.878 


0.0527 


0.0734 


0, 


.5010 





.9993 


1 


.801 


0, 


.8112 


0. 


,1438 


C12 


1 


.0024 


0, 


,5537 


31 


.901 


0.0497 


0.0713 


0, 


.5011 





.9961 


1 


.814 


0, 


.7593 


0. 


,1388 


C13 


1 


.0025 


0, 


,5537 


31 


.901 


0.0480 


0.0701 


0, 


.5010 





.9943 


1 


.819 


0, 


.7309 


0. 


,1360 


C14 


1 


.0022 


0, 


,5539 


31 


.934 


0.0465 


0.0689 


0, 


.5010 





.9933 


1 


.824 


0, 


,7048 


0. 


,1332 


C15 


1 


.0022 


0. 


,5537 


31 


.965 


0.0448 


0.0677 


0, 


.5011 





.9919 


1 


.829 


0, 


.6777 


0. 


,1303 


C16 


1, 


.0025 


0. 


,5533 


31, 


.971 


0.0431 


0.0664 


0, 


,5010 





.9905 


1 


.836 


0, 


,6507 


0. 


,1274 


C17 


1 


.0024 


0. 


,5533 


31 


.982 


0.0415 


0.0651 


0, 


,5010 





.9887 


1 


.841 


0, 


,6239 


0. 


,1245 


C18 


1 


.0021 


0. 


,5534 


32 


.020 


0.0399 


0.0639 


0, 


,5010 





.9873 


1 


.847 


0, 


.5983 


0. 


,1216 


C19 


1 


.0022 


0. 


,5532 


32 


.053 


0.0383 


0.0625 


0, 


,5010 


0, 


.9861 


1 


.852 


0, 


.5728 


0. 


,1187 


C20 


1 


.0025 


0. 


,5530 


32 


.057 


0.0366 


0.0612 


0, 


.5011 


0, 


.9850 


1 


.857 


0, 


.5488 


0. 


,1158 


C21 


1, 


.0024 


0. 


,5530 


32 


.072 


0.0351 


0.0598 


0, 


.5010 





.9834 


1 


.863 


0, 


.5242 


0. 


1129 


C22 


1, 


.0021 


0. 


,5525 


32 


.115 


0.0335 


0.0585 


0, 


,5011 





.9827 


1 


.867 


0, 


,5012 


0. 


,1100 


C23 


1 


.0021 


0. 


,5530 


32 


.152 


0.0320 


0.0572 


0, 


,5010 





.9815 


1 


.873 


0, 


.4788 


0. 


,1072 


C24 


1 


.0025 


0. 


,5528 


32 


.158 


0.0305 


0.0558 


0, 


,5011 





.9804 


1 


.877 


0, 


.4571 


0. 


,1044 



Table 9 — Continued 



Model Ml 


Ri 


i=i fi=i 

Pmax JKL 


ji=l 
spm 


Ma 


R2 


i=2 fi=2 
Fmax J RL 


ji=2 
spin 


(Mo) 


(10^ cm) 


(10^ cgs) 


(10^0 cgs) 


(Mo) 


(10^ cm) 


(10^ Cgs) 


(10^° cgs) 



C25 


1 


.0024 


0.5522 


32, 


.175 


0.0290 


0.0544 


0, 


,5010 


0.9793 


1 


.882 


0.4357 


0.1015 


C26 


1 


.0020 


0.5524 


32, 


.223 


0.0276 


0.0531 


0, 


,5011 


0.9786 


1 


.886 


0.4154 


0.0988 


C27 


1, 


.0021 


0.5523 


32 


.267 


0.0262 


0.0517 


0, 


,5011 


0.9775 


1 


.891 


0.3951 


0.0960 


C28 


1 


.0024 


0.5515 


32 


.281 


0.0248 


0.0503 


0, 


,5011 


0.9761 


1 


.895 


0.3755 


0.0933 


C29 


1 


.0024 


0.5515 


32 


.291 


0.0235 


0.0489 


0, 


.5011 


0.9753 


1 


.899 


0.3572 


0.0906 


C730 


1 


.0020 


0.5518 


32 


.342 


0.0223 


0.0476 


0, 


,5011 


0.9746 


1 


.903 


0.3393 


0.0879 


C731 


1 


.0020 


0.5516 


32 


.397 


0.0211 


0.0463 


0, 


,5011 


0.9739 


1 


.907 


0.3219 


0.0853 


C32 


1 


.0022 


0.5510 


32 


.425 


0.0199 


0.0450 


0, 


,5011 


0.9727 


1 


.911 


0.3049 


0.0827 


C33 


1 


.0029 


0.5505 


32 


.425 


0.0197 


0.0448 


0, 


,5011 


0.9727 


1 


.911 


0.3041 


0.0825 


C734 


1, 


.0021 


0.5511 


32 


.462 


0.0187 


0.0435 


0, 


,5011 


0.9722 


1 


.915 


0.2883 


0.0800 


C735 


1 


.0019 


0.5508 


32 


.522 


0.0176 


0.0422 


0, 


,5011 


0.9720 


1 


.919 


0.2735 


0.0775 


C736 


1 


.0018 


0.5508 


32 


.522 


0.0166 


0.0409 


0, 


,5009 


0.9711 


1 


.919 


0.2582 


0.0750 
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Table 10. Semi-detached DWD Sequence 'D'; Mtot = l.SM© 



Mtot 


q 


a 





Jtot 








VE 




(Ma) 




cm) 


(10-2 S-1) 


(10^° cgs) 


(IqSO gj.g) 










1.5048 


1.0000 


2.1010 


14.8472 


5.3866 


-1.8636 


1, 


.5 


X 10" 


-4 


1.5034 


0.9508 


2.1638 


14.1787 


5.4116 


-1.8566 


1, 


.8 


X 10" 


-4 


1.5045 


0.8994 


2.2336 


13.4444 


5.4491 


-1.8611 


9 


.5 


X 10" 


-5 


1.5034 


0.8507 


2.3102 


12.8287 


5.4797 


-1.8628 


1, 


.4 


X 10" 


-4 


1.5033 


0.8007 


2.3961 


12.1414 


5.5141 


-1.8736 


2, 


.2 


X 10" 


-4 


1.5034 


0.7504 


2.4919 


11.4431 


5.5460 


-1.8925 


3 


.1 


X 10" 


-4 


1.5034 


0.7004 


2.5987 


10.7442 


5.5702 


-1.9194 


4, 


.0 


X 10" 


-4 


1.5034 


0.6504 


2.7188 


10.0382 


5.5871 


-1.9572 


5, 


.5 


X 10" 


-4 


1.5032 


0.6004 


2.8550 


9.3217 


5.5914 


-2.0072 


7, 


.2 


X 10" 


-4 


1.5033 


0.5503 


3.0109 


9.0670 


5.5804 


-2.0741 


9 


.1 


X 10" 


-4 


1.5031 


0.5004 


3.1914 


7.8907 


5.5484 


-2.1587 


1, 


.2 


X 10" 


-3 



Table 11. Semi-detached DWD Sequence 'E'; Mtot = 1-0M( 



Mtot 
(Mo) 


Q 


a 

(10^ cm) 


n 

(10-2 S-1) 


Jtot 

(10^° cgs) 


-E'tot 
(IqSO grg) 






VE 




1.0030 


1.0000 


2.7109 


8.2477 


3.3579 


-0.6657 


1, 


.7 


X 10" 


-4 


1.0024 


0.9504 


2.7778 


7.9369 


3.3687 


-0.6642 


1, 


.6 


X 10" 


-4 


1.0019 


0.9001 


2.8641 


7.5710 


3.3850 


-0.6636 


2, 


.1 


X 10" 


-4 


1.0022 


0.8508 


2.9380 


7.2829 


3.3983 


-0.6659 


2, 


.2 


X 10" 


-4 


1.0024 


0.8003 


3.0161 


6.9993 


3.4073 


-0.6698 


2, 


.2 


X 10" 


-4 


1.0024 


0.7503 


3.1165 


6.6609 


3.4163 


-0.6751 


2, 


.6 


X 10" 


-4 


1.0024 


0.7003 


3.2294 


6.3129 


3.4212 


-0.6828 


3 


.3 


X 10" 


~4 


1.0023 


0.6502 


3.3569 


5.9551 


3.4204 


-0.6934 


3 


.9 


X 10" 


~4 


1.0023 


0.6003 


3.5026 


5.5862 


3.4123 


-0.7074 


4, 


.7 


X 10" 


-4 


1.0022 


0.5503 


3.6696 


5.2078 


3.3943 


-0.7254 


5, 


.6 


X 10" 


-4 


1.0021 


0.5003 


3.8649 


4.8169 


3.3644 


-0.7484 


6 


.8 


X 10" 


-4 


1.0020 


0.4504 


4.0955 


4.4134 


3.3184 


-0.7773 


8, 


.0 


X 10" 


-4 
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Fig. 1. — The mass-radius relationship is shown for spherical stars with our adopted ZTWD 
equation of state. Diamonds represent results derived using our three-dimensional SCF scheme 
applied to nonrotating, isolated configurations (see Tab l e [2l) : a sterisks show previously published 
results for the same equation of state taken f romlHachisul ( 1986al ): the solid curve shows the analytic 
mass-radius relation, Eq. ()A14p . derived by Nauenbe^ ( 19721 ). 
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Table Al. Physical Constants 



Constants'^ 


This Paper*^ 


Chandrasekhar fl967l^ 


(1) 


(2) 


(3) 


c (cm s~^) 


2.99792 X 10+^° 


2.9978 X 10^° 


h (erg • s) 


6.62608 X 10^27 


6.62 X 10-27 


me (g) 


9.10939 X 10^28 


9.105 X 10"28 


rrip (g) 


1.67262 X 10^24 


1.672 X 10-24 


(g) 


1.66054 X 10^24 




G (cm^ g~^ s~^) 


6.6726 X 10-8 


6.62 X 10-8 


Mo (g) 


1.9891 X 10^3 


1.985 X 10^3 


Rq (cm) 


6.955 X 10^° 


6.951 X 10^° 


A (dynes cm~^) 


6.00228 X 10^2 


6.01 X 10^2 


B^~^ (g cm~^) 


9.81011 X 10^ 


9.82 X 10^ 


^i/ie (cm) 


7.71395 X 10^ 


7.705 X 10^ 



'^Speed of light, c; Planck's constant, h; mass of the electron, 
TTie; mass of the proton nip] atomic mass unit, m^; universal 
gravitational constant, G; solar mass, Mq\ solar radius, Rq; as 
used in the ZTWD equation of state ([I]), j4 = irm^c' /3h^ and 
Bii~^ = 8-Km^c^mp/3h^; the characteristic WD length scale, 

hl^e = {2A/7TGy/^{fle/B). 

^Drawn fromlCoJ kood ). 



^Drawn from Appendix I, Table 32 of lChandrasekhaij (|l967l ). 
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Fig. 2. — Schematic diagram illustrating the equatorial-plane structure of a binary star system. 
The primary star, on the left, has a mass Mi, a radius and a central density pj^xi the secondary 
star, on the right, has a mass M2 < Mi, a radius R2, and a central density p^ax- The centers of 
mass of the two stars (points labeled Oi and O2) are separated by a distance a = wi + W2, and 
their distances from the center of mass of the system are, respectively, wi and ■W2. The points 
labeled Oa and Op identify, respectively, the outer edge and inner edge of the secondary star. 
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Fig. 3. — Density contours in the equatorial plane are shown for eight separate ZTWD binary models with increasing separation 
along inspiral sequence 'A' (Mtot = 1-5Mq; q = 1). Labels in the upper-right-hand corner of each panel identify each model by 
number according to its corresponding position along the sequence as itemized in Tables H] and [5l 
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Fig. 4. — Density contours in the equatorial plane are shown for eight separate ZTWD binary models with increasing separation 
along inspiral sequence 'B' (Mtot = 1-5Mq; q = 2/3). Labels in the upper-right-hand corner of each panel identify each model 
by number according to its corresponding position along the sequence as itemized in Tables [6] and [71 
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Fig. 5. — Density contours in the equatorial plane are shown for eight separate ZTWD binary models with increasing separation 
along inspiral sequence 'C (Mtot = 1.5Mq; q = 1/2). Labels in the upper-right-hand corner of each panel identify each model 
by number according to its corresponding position along the sequence as itemized in Tables [8] and [H 
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Fig. 6. — (Top panel) Normalized total angular momentum, Jtot /'^norm, {second panel) normalized 
total energy, StotZ-E'normi {third panel) the Roche-lobe filling factor, /rl, for the secondary (aster- 
isks) and primary (diamonds) stars, and {bottom panel) the normalized orbital angular velocity, 
O/f^K) are plotted as a function of binary separation for models Al through A35 along inspiral 
sequence 'A' (Mtot = 1-5Mq; q = 1). Data for the individual models is drawn from Tables H] and [Sj 
the separation a is labeled in units of 10^ cm along the bottom axis and as a ratio to i?o.75 along the 
top axis. The solid curves in the top two panels display the analytic functions Jpm{(i) / Jnorm and 
-f'pm (a) /-Bnorm given, respectively, by Eqs. ([lOj) and for a "point-mass" sequence of the specified 
total mass and mass ratio; and the dot-dashed curves display the analytic functions Jss(a)/«/norm 
and E^Ja) I Er^nrrr, aoDroDriate for a "sDinning spheres" sequence given, resoectivelv, bv Eos. (I4T]) 
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Fig. 7. — Same as Fig. [6] but for models Bl through B49 afong the inspiral sequence 'B' (Mtot = 
I.5M0; g = 2/3), as tabulated in Tables [6] and O along the top axis, the separation a is labeled as 
a ratio to Ro.eo- 
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Fig. 8. — Same as Fig. [6] but for models CI through C36 along the inspiral sequence 'C (Mtot = 
1.5Mq; q = 1/2), as tabulated in Tables [8] and [9l along the top axis, the separation a is labeled as 
a ratio to -Rq.so- 
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Fig. 9. — System parameters a, Q, Jtot, and i?tot at contact as a function of mass ratio, q, for 
DWD systems having a total mass of 1.5Mq. Solid curves in the top three panels sho'w predicted 
behavior based on Nauenberg's (1972) and Eggleton's (1983) approximate, analytic expressions as 
discussed in the text. 



